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Abstract 


The robustness of detection filters applied to the detection of actuator failures on a 
free-free beam is analyzed. This analysis is based on computer simulation tests of 
the detection filters in the presence of different types of model mismatch, and on 
frequency response functions of the transfers corresponding to the model mismatch. 

The robustness of detection filters based on a model of the beam containing a large 
number of structural modes varied dramatically with the placement of some of the 
filter Doles. The dynamics of these filters were very hard to analyze. 

The design of detection filters with a number of modes equal to the number of 
sensors was trivial. They can be configured to detect any number of actuator 
failure events. The dynamics of these filters were very easy to analyze and their 
robustness properties were much improved. A change of the output transformation 
allowed the filter to perform satisfactorily with realistic levels of model mismatch. 
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Chapter 1 
Introduction 

It is expected, with the advent of the space shuttle, that large space 
structures serving a variety of mission requirements will be assembled in space in 
the near future. Some of the structures being contemplated include antennas, 
reflectors and solar power satellites. They are characterized by their light weight 
and very large sizes, which in turn result in low frequency bending properties with 
very little structural damping. 

The active control of large space structures will include station keeping 
control, attitude control, shape control and vibrational damping. The component 
unreliability issue is most important in the last of these applications where 
hundreds of sensors and actuators will be needed in order to actively damp the 
many vibrational modes to assure mission success. The large number of 
components in the system result in a high probability of a component failure in a 
short period of time, even assuming components having very optimistic mean time 
between failures For example, a system with 200 components, each having a time 
between failures of 100,000 hours, is expected to have 17 component failures a 
year (7). Even if these control systems are serviced in orbit, it is desired that its 
service interval be as long as possible, at least one year. Therefore, it will be 
necessary to design control systems that can tolerate the number of component 
failures occurring during the service interval. This, in turn, will require systems 
that can detect and identify a component failure, and will require control systems 
that can be reconfigured to perform without the failed component in some optimal 
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fashion. 

The failure detection and identification systems (FDI) depend on hardware 
redundancy or analytic redundancy in order to detect and identify a failed 
component. An example of hardware redundancy to detect and identify sensor 
failures is the one that involves triplication of sensors, where the failure of a 
component is detected when there is a discrepancy between the signals of two 
sensors and comparison with the third determines which of the two has failed. 
Because of the large weight and cost resulting from hardware redundancy, these 
systems are not attractive for applications in large space structures. The FDI 
systems based on analytic redundancy use knowledge of the plant dynamics to 
make the detection and identification of the failure. Some of these systems require 
specification of the failure modes ahead of time. 

The Failure Detection Filter is a closed loop method, based on analytic 
redundancy, to detect and identify actuator and sensor failures, and changes in 
plant dynamics. A very attractive property of the detection filter is that it does 
not require specificaton of the mode of failure. 

The Failure Detection Filter is a linear filter having the same configuration as 
an observer but with an additional constraint in the design of the gain matrix 
D. The output error, the difference between the system measurements and the 
filter estimate of those measurements, is zero when the system is operating 
nominally (assuming no model mismatch). When a failure occurs, the output error 
becomes nonzero, indicating that a component has failed. Moreover, the additional 
constraint on the gain matrix D keeps the output error in a single direction in the 
output space. The failed component can then be identified by this direction. 

The failure detection filter was first proposed by Beard [lj for deterministic 
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systems. Jones [5] explained detection filter theory using a more geometric 
approach and expanded it to stochastic and sampled data systems. Detection 
Filter Theory was applied by Carignan [2] in the context of flexible space 
structures, where he considered actuator and sensor failure events, and analyzed 
the performance of the filter in the presence of unmodeled modes. He developed 
the computer software used in this thesis to design detection filters and to simulate 
the dynamics of the plant with the filter. 

The goal in this thesis, is to improve the performance of detection filters in 
detecting and identifying actuator failures in an undamped flexible structure, in the 
presence of model mismatch. To achieve this goal, the filter performance will be 
evaluated in the presence of unmodeled modes and different levels of parameter 
errors. In addition, the transfer properties of the filter will be analyzed and the 
results compared with the simulations. 

In Chapter 2, the main results of detection filter theory are introduced along 
with analytical procedures for the filter design. Chapter 3 describes the model of 
the undamped flexible beam, shows some properties of the detection filter design 
for actuator failures, and simulates the filter with no model mismatch. 

In Chapter 4, the detection filter performance in the presence of model 
mismatch is analyzed using simulations and frequency response functions. Ways of 
improving the filter performance derived from this analysis are presented in the 
chapter. Chapter 5 contains the conclusions and recommendations for further 
research. 
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Chapter 2 

Failure Detection Filter Theory 

The Failure Detection Filter is a. linear filter having the same 
configuration as a full order observer (Fig. 2-1). It includes a model of the nominal 
system and is driven by the same actuator commands as the actual system. The 
output error, the difference between the system measurements and the filter 
estimate of those measurements, is multiplied by the gain matrix D and fed-back to 
the filter. 

In a full order observer, the only constraint placed upon the gain matrix D is 
that the eigenvalues of the matrix (A - DC) must have negative real parts. When 
this is the case, the observer is stable and the state and output errors become zero 
(except for noise and other unmodeled effects) during nominal operation. When a 
component failure occurs, the output error will become nonzero. Therefore, any 
stable observer can detect when the system has failed. What distinguishes a F ailure 
Detection Filter from such an observer is that the gain matrix D has the additional 
constraint of restricting the output error due to a particular failure to a single line 
in the output space. The direction of that line is used to identify which component 
has failed. 

Pole placement in Detection Filter design can be used not only to assure the 
stability of the filter, but to improve the performance of the filter in the presence 
of unmodeled system dynamics. Therefore, complete control of the Detection Filter 
poles is required. This requirement together with the output directionality 
constraint are the bases of Detection Filter Theory and are present in the definition 
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Figure 2-1: The Failure Detection Filter 
of Failure Detectability. 

Detection Filter Theory consists mainly of two parts: 1) a set of necessary 
and sufficient conditions that must be satisfied for the solution of a particular 
detection filter problem to exist, and 2) a set of algorithms to achieve the design. 

In this chapter, the state and output error equations are presented for both 
nominal and failed system operation. The concept of Failure Detectability will 
then be introduced, and important results concerning necessary and sufficient 
conditions will be stated. Finally, an algorithm for filter design based on 
annhilating minimal polynomials will be presented. 
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2.1 The Failure Detection Filter Problem 

The model that describes the plant dynamics must be linear, time-invariant 
and observable. It is not necessary that the system be observable but only that the 
model be observable. An observable model can be obtained by removing the 
unobservable dynamics from it. This does not reduce the effectiveness of the FDI 
system because failures that affect only the unobservable part of the plant 
dynamics cannot be detected. 

The unfailed system is represented by the equations: 

i(<) = Ax(l) + B»(!) (2.1) 

y(f) = C x(t) 

where 

x(f) is the n dimensional state vector, 
u(f) is the m dimensional control vector, 
y(t) is the p dimensional measurement vector. 

The structure of the Failure Detection Filter is shown in (2-1). Due to the 
similarity between detection filters and observers, the filter state vector x is called 
the state estimate and f is called the measurement estimate. 

The filter equations are: 

W) = A x(t) + B u c (/) + D [ y(l ) - y(<) ] (2.2) 

y«) = C x(l) 


The state error is defined as: 
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e ( f ) = x (<)-*(0 

The measurement error is defined as: 

<(<) = r( ')-*<<) 

The state and measurement equations are obtained by subtracting equations 
(2.1) from (2.2). In the normal mode, u(f) = u c (f)> and 

£(<) — (A - DC) e(f) (2.3) 

e(t) = C e(/) 

If the eigenvalues of the matrix (A - DC) have negative real parts, the state 
and measurement errors will become zero in the steady state. This is true, 
however, only when the system is operating nominally. When a component failure 
occurs, the state and output error equations (2.3) are no longer valid since 
equations (2.1) no longer represent the true system. The new state and output 
error equations will depend on the type of component that failed. 

Consider the failure of actuator ;. The true control u (t) that is being applied 
to the plant is 

«(0 = «<.(•) + « mj lit) 

where u c ^) is the desired control, e m j is a unit vector in the jth direction of 
dimension m, and i\(t) is a scalar time function equal to the difference between the 
commanded input to actuator j and the actual control being applied by this 
actuator. 


The new system equations are 


*(<) = A x(/) + B u c (Z) + b y ri(t) 

7(0 = C x(/) 
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(2.4) 


where bj is the column of B corresponding to the jth actuator. 

The new state and output error equations with actuator j failure are obtained 
by combining (2.2) and (2.4) 

i(f) = (A- DC) e(<) + !,;,)(() (2.5) 

<(<) = C *(() 

The goal of Detection Filter design, in this example, is to find a gain matrix 
D such that e(f) maintains a fixed direction in the output space when actuator j 
fails (and at the same time assign almost arbitrarily the eigenvalues of matrix [A- 
DC]). Since only directionality is important in the identification of the failed 
component, the vector b >. is the only element in equation (2.5) that characterizes 
the failure of actuator j. For this reason, the vector b^ is called the event vector 
associated with actuator j 'ailure. The forcing function rj(t) depends on the mode 
of failure of the actuator, and knowledge of this signal (or of the failure mode) is 
not needed in the use of a detection filter. This characteristic makes Failure 
Detection Filters very attractive for failure detection and identification. 

For the rest of the chapter, the detection problem is generalized to the 
detection and identification of the component failure associated with the event 
vector f, and whose corresponding state and output error equations are 

6(0 = (A-DC)e(f) + ft?(0 (2.6) 

€(/) = e(<) 

The discussion in this section, so far, can be summarized by the following 



definition of detectability found in Beard [1]. 

Definition: The failure event associated with event vector f in (2.6) is 

detectable if there exists a matrix D such that 

(1) C ejl) maintains a fixed direction in the output 
space (where e s (f) is the settled-out solution 

of (2.6) ) and, 

(2) at the same time, all eigenvalues of (A - DC) can be 
specified almost arbitrarily. 

Condition (1) is the distinguishing feature of the failure detection filter and 
serves to identify which component has failed. Condition (2) is needed to assure 
that the matrix (A - DC) is stable so that during nominal operation, equation (2.3), 
the initial state and output errors die out to zero, and in the presence of component 
failure, equation (2.6), the output error is unidirectional. Condition (2) also allows 
the designer of a failure detection filter to adjust the bandwidth of the filter to 
suppress the effects of unmodeled system dynamics, as will be shown in Chapter 5. 

Given a failure event, the detection filter problem consists of: determining 
whether the event is detectable and, if it is detectable, to find the gain matrix 
that makes the output error unidirectional and assigns the desired poles to the 
filter. 

In the next section, the detection filter problem is solved for both fully 
measurable systems (rank C = n) and partially measurable systems (rank C < n). 
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2.2 Fully Measurable Systems 

A fully measurable system is a system whose state vector x(f) can be 
solved uniquely given the measurement vector y(t), for any time t. Therefore, a 
system is fully measurable if and only if 

rank C = n 

this implies that the number of sensors p > n. 

The solution of the detection filter problem for this type of system is very 
simple. Choose a gain matrix D such that 

(A - DC) = - c I (2.7) 


where I is the identity matrix and a is a positive scalar constant. Then all the 
eigenvalues of (A - DC) are -<r, and the resulting filter is stable. The matrix D 
that satisfies (2.7) is given uniquely by 

D = (A + <r I) C' 1 (2.8) 

if rank C = n = p, and non-uniquely by 

D = (A + or I) (C 1 ^)- 1 C T (2.9) 

if rank C = n < p. 

To prove that this choice of D satisfies condition (1) of detectability, a result 
trom linear systems theory will be used. When the failure associated with event 
vector f occurs, the state error in (2.6) will be driven by f t j(t) within a subspace 
of R n called the controllable space of f. This subspace is given by the range space 
of the matrix W f defined as 
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W f = 


f (A-DC)f . . . (A -DC)"-^ 


( 2 . 10 ) 


Since £ = C e, the output error will lie in a subspace of given by the range 
space of the matrix CW f . 

For thi3 choice of matrix D, the matrix CW f becomes 


CW f = 


Cf -<rCf <r 2 Cf . . . 


and the range space of CW f is simply the one-dimensional space spanned by the 
vector Cf. Therefore, when the failure occurs, the output error will remain in a 
single direction given by the vector Cf. Hence, the gain matrix D given by (2.8) or 
(2.9) also satisfies condition (1) of detectability. Any matrix D which satisfies this 
condition is called a detector gain for f. 

From the preceding proof, note that this choice of D (2.8) or (2.9) is a 
detector gain for any event vector f in R”. This implies that a single detection 
filter can be used to detect and identify every failure associated with every event 
vector f in R” (except that it wouldn’t be able to identify two failures that have the 
same event vector f). 

The solution of the failure detection filter problem in the case of fully 
measurable systems is trivial. Choosing (A - DC) = -or I assured the stability of 
the filter and also produced unidirectional output errors in the event of component 
failures. 
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2.3 Partially Measurable Systems 

In a partially measurable system rk C < n and therefore, the state vector 
x(<) cannot be solved uniquely given the measurements y(f). In the previous section 
it was shown that when a system is fully measurable, a single detection filter can 
be designed to detect and identify every event f in R n . This property is lost when 
the system is partially measurable, and more than one detection filter may be 
needed to detect a given set of failure events. Moreover, in such detection filters, 
there are subspaces of R" containing event vectors that produce the same 
unidirectional output. These subspaces are called Detection Spaces and are very 
important in the analysis and design of detection filters, for both single and 
multiple detection. 

2.3.1 Single Failure Detection 

This section explores the detection problem for a single failure event. The 
definition of detection space is introduced and is later used to prove the 
detection theorem. Algorithms for the design of detection filters, based on the 
proof of the detection theorem, are then discussed. 

2.3.1. 1 The Detection Space 

Assuming that (A,C) is an observable pair, Jones defines the detection space 
for the event f in the following way: 

Definition 2.3-1 Assume that Cf is not zero. The detection space for f is 
denoted by R f and is the direct sum: 

i) Rf = f 0 

where R f C R” is the largest subspace which satisfies the two 
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conditious: 

ii) R f C rf(C) 

iii) A R f C Rf 

This definition is motivated by condition (1) of detectability, and serves to 
derive a formula to generate detector gains. By condition ii) and iii) of Definition 
2.3-1: 

(A - DC) R f C A R f C R, 

(A - DC) R f C R f (2.11) 

Assuming that Cf 0, condition i) and ii) imply that: 

C R f = Cf 0 C R f = Cf 0 0 = Cf 

dim (C R f ) = dim (Cf) = 1 (2. 12) 

Now, assume that D is chosen such that: 

(A - DC) f = £ € Rj (2.13) 

for some arbitrary £ in Rf . From (2.11), (2.13) and i), it follows that R f is an 
invariant space with respect to (A - DC) for this choice of D since: 

(A - DC) R f = (A - DC) (f 0 R f ) C R f 

(A - DC) Rf C Rf (2.14) 

Since Rf contains the event vector f, and since Rf is invariant with respect to 
(A - DC) (for this choice of D), the controllable space of f, spanned by the columns 
of the matrix defined in (2.10), is a subspace of Rfi 
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(2.15) 


C t C Rf 

where C f is the controllable space of f. 

Then, from (2.12) and (2.15): 

CC r =Cf (2.16) 

dim (CC f ) = 1 

and the failure associated with f generates unidirectional output errors along Cf. 
Therefore, the gain matrix D chosen in (2.13) is a detector gain. Moreover, by 
the same argument that led to (2.16) and since (A • DC,C) is an observable pair, 
any event vector in produces the same unidirectional output Cf. Jones [5] calls 
this property detection equivalence and proved that the detection space R f 
contains all the event vectors that are detection equivalent to f. 

Since, as it was proved, the gain D satisfying (2.13) is a detector gain, this 
equation can be solved for D to obtain an expression for the detector gain of f. 

Equation (2.13) can be written as: 

D Cf = A f - £ (2.17) 


(2.18) 


where £ e Rp and its solution is 

D = D p + D h 

where 

Dp = (a r- Merger]- 1 (cr) T 

D„ = D' 


I - Cf [ (Cf) T Cf]’ 1 (Cf) T 
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with £ e E f and D’ an arbitrary matrix with the same dimension as D. 

The detector gain of f is given by (2.18) with any ( c Rj and arbitrary D' 
The eigenvalues of (A - DC) will depend on the choice if £ *■ Rj and D’, and the 
question is whether all the eigenvalues of (A - DC) can be arbitrarily specified by 
the proper choice of £ £ Rj and D\ This question is addressed in the next section. 

Up to this point, it has been assumed that Cf ^ 0 . If Cf = 0 , it can be 
shown that the definition of detection space and all the results derived so far, can 
be extended by replacing f by A“ f, where u > 0 is the smallest integer, such that 
CA U f 0 . To simplify the analysis, it will be assumed for the rest of this chapter 
that Cf 5 ^ 0. 

2.3. 1.2 Detection Theorem 

In this section, it will be shown that the freedom represented by the choice of 
£ e R f and D’ in equation (2.18) is enough to assign all the eigenvalues of (A - DC) 
arbitrarily. 

First, a theorem proved by Beard and Jones is stated. 

Theorem 2.3.1 Let dim (Rf) = v { . There exists a unique vector g < R f , 
called the detection generator of Rp such that: 

i) A k g e R f k = 0 , 1 , 2 , . . . , (t> f - 2 ) 

ii) C A v f l g = Cf 

This theorem implies that the vectors g, A g, . . . , A u f ^ g are a basis for R f . 
Since f e Rp f can be expressed as 

f=Q,g + 0i 2 Ag + ... + a , A v (- J g + AV 1 g 


(2.19) 
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where *= 1 by ii). 

Also, since ( e R p ( can be written as 


( = a 1 g + a 2 A g + . . . + a^ AV 1 g 
Write (2.13) as 


D Cf = A f - £ 


and use (2.10) and (2.20) to get 

A f - ( = - a, g + (a, - » 2 ) A * + (a 2 - » 3 ) A 2 g + . . . 

■ + K,-l - a Uf ) A Vf l * + ^'8 

Define = -a r p 2 = Qj - a 2 , . . . , p = a j - a^, and write (2.22) as 
A f *£ = P! g + p 2 A g + . . . •}' p Uf AV 1 g + A u f g 
Equation (2.21) can now be written as 

DCr=p 1 g + ... + p l)( AV'g + A t ’tg 
It is easy to prove, using Theorem 2.3.1, that 
A k g = (A - DC) k g for 0 < k < i/ f 
Therefore, (2.24) can be expressed as 

0 * Pj S + P 2 (A - DC) g + . . . + p Vf (A - DC)V l g + A v f g - D Cf 
and, since Cf = C A u f '* g, using (2.25) we get 

0 = P! 8 + P 2 ( A - DC) g + . . . + p Vf (A - DC) v f * ! g +. (A - DC) u f g 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


(2.24) 


(2.25) 


(2.26) 
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Write (2.28) as 

V>(A - DC) g = 0 

where the polynomial ^-) is defined 

Mi) = ?i + p 2 7 + • • • + p Uf n v t ‘ l + i V{ 

The polynomial ^( ) is the minimal annihilating polynomial of g with 
respect to (A - DC) (5). Therefore, the eigenvalues of (A - DC) associated with the 
detection space R f are given by the roots of 

m = o 


or 


Pj + P 2 X + . . . + p Uf \ v r l + \ v t = 0 (2.27) 

For each set of eigenvalues of (A - DC) associated with R f , th°re is a unique 
set of coefficients pj, p 2 , . . . , p a unique set of coefficients a^ a„, . . . , a^ and, 
therefore, there is a unique £ € R^ Hence, the detection space eigenvalues of (A - 
DC) can be assigned arbitrarily by the proper choice of ( t R f in (2.18). Condition 
ii) of detectability will be satisfied only if the remaining n - Uj eigenvalues can be 
assigned arbitrarily by choice of D’ in (2.18). 

By replacing D from (2.18) into (A-DC) get 


A - DC * A’ - D C’ 


where 


A’ = A - DC 

p 


and 
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C’ = 


I-Cf [(Cf) T Cf]- 1 (Cf) T 


(2.28) 


Matrix D from (2.18) makes R f invariant with respect to (A - DC) and (A’- 
D’C’). Since D’ is an arbitrary matrix, Rj is also invariant with respect to A’: 


A’ Rf C Rj 


(2.29) 


Since: 


C' R f = C’ (f © R f ) = C’f © C’R f = 0 
the detection space R f is in the null space of G”: 

R f C *?(C ) (2.30) 

From (2.29) and (2.30), it can be concluded that the detection space R f is a 
subspace of the unobservable space of the (A’, C’) pair: 

R f C »?(M) (2.31) 

where M is defined as 

M = 

Since: 


C’ 

C’A’ 

CA’”* 1 


C’ ri( M) = 0 


using equation (2.28) get 



C ri( M) = Cf 
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(2.32) 


But since »f(M) is an invariant space with respect to A’ and to (A - DC), 
equation (2.32) implies that 

l(M) C R, (2 33) 

Hence, by (2.31) and (2.33): 

R, = ij(M) (2.34) 

Given an unobservable pair (A’, O’), the theory of linear systems says that all 
of the eigenvalues of (A’ - D’C’) can be placed arbitrarily by choice of D’ except 
those eigenvalues associated with the unobservable space of (A 1 , C’), which are 
fixed. But since the unobservable space of (A’, C’) is the detection space of R f , the 
fixed eigenvalues of (A’ - D’C’) are the ty eigenvalues associated with R f , and were 
assigned by the choice of £ e R f in (2.18). Therefore, the matrix D’ in (2.18)) 
assigns arbitrarily the remaining n - eigenvalues of (A - DC) and condition ii) of 
detectability is satisfied. 

The results of this section are summarized with the statement of the 

Detectability Theorem. 

Theorem 2.3.2 Every vector in the state space R” is detectable if and only 
if (A,C) is observable. 

2.3. 1.3 Algorithms 

The algorithms presented in this section are based on the proof of the 
detectability theorem, and were suggested by Beard (1). 

(1) Detection Space Algorithm 
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This algorithm is used to obtain the detection space R f and its dimension 
of the event f. From (2.34), the detection space R f is in the null space of the 
matrix M defined in (2.31). The matrices C’ and A’ present in this definition were 
introduced in (2.28), where C’ depends only on f and A’ is a function of D . The 
detector gain D p is defined in (2.18) for any £ t R f . Therefore, the detection space 
R f is in the null space of any M that results from any choice of £ c R f in (2.18). 
Since, from definition(2.3-l), the only vectors £ known to be in R f are the ones in 
the subspace spanned by f, the vector £ = 0 e R f is chosen and the A’ resulting 
from this choice is called K f : 

K f = A - A f [(Cf) T Cf]' 1 (Cf) T C (2.35) 

Hence the algorithm to obtain R f and is 

R f = *?(M f ) 

= n - rank M f 

where 



C’ 

C’K. 

• * 


CKf 1 


(2.36) 


with C’ given by (2.28) and K f by (2.35). 

(2) Detection Gain Algorithms 

These algorithms perform the actual detection filter design. They calculate 
the gain matrix D that satisfies condition i) of detectability, and assigns the matrix 
(A - DC) the desired eigenvalues. 
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The gain D given by equation (2.18) is a detector gain, and therefore satisfies 
condition i) of detectability. The choice of £ e R f assigns the eigenvalues of (A - 
DC) associated with R f , and the choice of D’ assigns the remaining n - t f 
eigenvalues associated with the completion of the state space R n . 

(a) D p Algorithm 

This algorithm calculates D p in (2.18) given a specified set of v f eigenvalues to 
be assigned to (A - DC), and that will be associated with R f . If the eigenvalues 
are specified as the roots of the polynomial 

s v t -I- p„ s v t' l + . . . + p.> s + p, = 0 

then, as it was proved, the matrix (A f - 0 must satisfy equation (2.23). Hence, by 
substituting (2.23) into (2.18) get 

D p = q [(Cf) T Cf]' 1 (Cf) T (2.37) 

where 

q = Pi g + P 2 A 8 + • • ■ + P Uf A ” rl 8 + AVf 8 

The detection generator g can be found by first using algorithm 1) to obtain 
the detection space R f and then apply the property of g stated in Theorem 2.3.1 to 
get g ( R f . A convenient way of executing these two steps is using orthogonal 
reduction, as suggested by Beard (lj. 

Once g is obtained, the matrix D p is simply calculated using (2.37) and is 
unique for a given set of eigenvalues. 

(b) D h Algorithm 

This algorithm assigns the remaining n - eigenvalues of (A - DC) by 
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choosing the appropriate D’ in (2.18) to calculate D^. These n - v f eigenvalues are 
associated with the completion of the state space R n , and their specifications do 
not determine a unique D’ (and a unique Dpj and D). Therefore, given an event f 
and a complete set of eigenvalues to be assigned to (A - DC), the solution of the 
detection problem will not be unique. For simplicity, the algorithm presented here 
is based on the same procedure used to assign the eigenvalues to R f and obtain D p . 
For this reason, this algorithm adds certain constraints to the detection problem 
and generates a unique D H that together with D p give a unique gain D as the 
solution of the detector problem. 

The matrix D H is given by (2.18), where D’ can be considered as the gain 
matrix of the system (A’, C’): 

A - DC = A’ - D’G" 
with A’ and C’ given by (2.28). 

Beard showed that there are a set of vectors w ; that have the same property 
with respect to (A’, C’) as the detection generator g with respect to (A, C): 

C'A’ P Wj=0 for 0 < p < q . -1 

C’A’V 1 w 0 

There are (rank C - # of events) of these vectors. They and the 
corresponding q^ can be obtained as a byproduct, when applying orthogonal 
reduction to find the null space of M in (2.31) or in (2.36) (Beard showed that M 
— M f ). 

The vectors Wy, A’wy, . . . , A ,q ; ~ l span a subspace called completion 
space Cy These spaces are the equivalent of R f to (A’, C’). It can be shown that 
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the spaces R f , Cj , C 2 , are all nonintersecting and 

R" = Rf 0 Cj e C 2 e . . . 0 C ( (2.38) 

where 

C t - H Cj — 0 , * j 

Rf n Q = o 

and 

l 

T,'ij =n - v t 

j-i 

This algorithm assigns q. eigenvalues to (A’ - D’C’) by making Cj an 
invariant space with respect to (A’ - D’C’). Therefore, the same algorithm used to 
find D p is used to get D’. Equation (2.37) is modified by substituting A’ for A, w. 
for g and C’A’ q i' 1 for Hence, the D’ that makes all C - invariant and assigns 
the remaining n - eigenvalues of (A’-D’C’), is given by 

D’ = * j(C’ W) T C’ W]* 1 (C’ W) T (2.39) 

where 

* = ’ • • * ’ 

= p{ w ; . + . . .+p j A’V 1 w y + A’V w ; . 

W = A’^l* 1 Wj , . . . , A’V 1 w, 

Beard showed that W can also be calculated as 
W= Kjl - 1 Wj , . . . , KJ/- 1 w, 
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where K f is defined in (2.35). 

2.3.2 Multiple Failure Detection 

This section deals with the problem of designing a single detection filter to 
detect a set of failure events {fj. . .f f }. This set of events is mutually detectable 
by a single detection filter, if there exists a D that satisfies the conditions of 
detectability for all the . For single failure detection, if (A,C) is an observable 
pair, any event f e R" is detectable. But for multiple failure detection, complete 
observability is not a sufficient condition. The Group Detection Theorem 
states a necessary and sufficient condition for a set of output separable vectors, 
as defined below, to be mutually detectable. 

Definition: The vectors {fj. . .f f } are output separable if 
rk CF = r 
where 

F = A U 1 fj , . . . , A u r f r (2.40) 

with tr > 0 defined as the smallest integer, such that CA“i f t 7^ 0. 

It can be shown that if the set of events are output separable 

Ef.nR r .= 0 , iVj (2.41) 

Before stating the Group Detection Theorem, the following definition is 
necessary. 

Definition: The group detection order of the set {fj . . .f r } is defined as 
the dimension of the null space of Mp, (n - rk Mp) where Mp is defined as M f in 
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(2.36) with f replaced by F in C’ and K f . 

Group Detection Theorem The output separable vectors {fj . . f r } are 
mutually detectable if and only if the sum of the individual detection orders of the 
set f t is equal to the group detection order. 

2.3.2. 1 Algorithms 

Given a set of event vectors {fj. . ,f r }, the algorithm must first determine 
whether the events are output separable by using (2.40). If they are output 
separable, it must then determine whether they are mutually detectable by 
applying the group detection theorem. 

If the events are output separable and mutually detectable, the 

procedure to calculate D is the same as the one described in the previous section for 
single failure detection, but with f replaced by F. 


D - D p + D h 


where 


D p = Q [(CF) t CF ]* 1 (CF) t 


(2.42) 


with 


Q = 


qt ’ • • > q r 


q t - = Pi g, + • • • + Py . AV 1 g| . + A v i gj 

where g^ and are the detection generator and detection order of the event T 


D h = D' 


I - CF [(CF) t CF ]- 1 (CF) t 


where 



-34- 


D’ = * [(C’ f W) t C’ F W j- 1 (C’ f W) t 

with 

W= K q l #1 w 1 ,...,K q l' 1 w / 

and 

V — i>f 

lp. = Pi W j + • ■ ■ + P^. A ' q i ' 1 W ;+ A ’^' W j 

where C” F and A’ are defined in (2.28), with A’ calculated using D p in (2.42) and 
C’ F calculated using F in place of f. 

Beard (1] shows a method to obtain the vectors w . and the corresponding q ., 
that consists of applying orthogonal reduction to find the null space of the matrix 

m f . 
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Chapter 3 

Dynamics of Flexible Space Structures 

Flexible structures are characterized by very low natural damping and by 
infinitely dimensioned plant dynamics. The design of detection filters for flexible 
space structures must, therefore, be based on a reduced order model of the system. 
Moreover, the parameters of the reduced order model are not exactly known, and 
the effect of parameter errors and of unmodeled dynamics can be very serious. 

In this thesis, the dynamics of a uniform flexible beam are used to analyze the 
performance of detection filters in the presence of model mismatch. The model of 
this flexible beam corresponds to an experimental beam that was assembled at 
NASA Langley Research Center. In this chapter, the model of the experimental 
beam will be described along with the design of detection filters for actuator failure 
events. The detection filter wiil then be evaluated with no model mismatch. 

3.1 The Flexible Beam Mode! 

The beam is made of aluminum (M = 0.502 slugs) and is twelve feet long 
with a 6” x 3/16” cross section [4j. It is equipped with four force actuators and 
four colocated displacement sensors (Fig. 3-1), that apply forces and measure 
displacement along the same direction in the vibration plane. 

The partial differential equation that describes the dynamics of the undamped 


beam is 
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12 3 4 



Moss =0.5031 slugs 


Figure 3-1: The Simulated Uniform Beam 


d 4 y d“y 

El a? +m d? = ^' l) !3 '" 


where £ = position along the beam 


y = displacement 


E = modulus 


m = mass per unit length 


I = cross section inertia 


f = distributed force per unit length 


The distributed force f(e , /) given by the four point actuators is 


f ( f . 0 = L 6 ( ( • e j ) U j0) 
;' = i 


(3.2) 



-37- 


where (.■ is the position of the jth actuator and is the control being applied by 
it. 

The solution of the partial differential equation is 
00 

y( f . 0 = H H*) V'/O (3 3i 

»— i 

where are the mode shapes and ^ t {t) are the modal amplitudes. The mode 
shapes <^e) comprise an orthogonal and normalized set of functions over (0, /), and 
the modal amplitudes ip t (t) satisfy an uncoupled set of differential equations 

1 * 

J— + “i H*) = - J2 » = 1, 2, . . . oo (3.4) 

;= i 

where M = m X 1 is the mass of the beam and u- is the natural frequency of the 
ilh mode. 

The orthogonal mode shapes <p i and the natural frequencies were obtained 
for the first ten modes at Langley, by performing a finite element analysis [4]. In 
Table (3-1) the natural frequencies and mode shape values at positions 0.5 ft., ‘2.5 
ft., 6.0 ft., and 9.5 ft., are given for the first eight modes - two rigid body modes 
(translation and rotation) and the first six bending modes. 

The kth sensor y k is located at and from equation (3.3) 

00 

y k (<) = E <*.K) H>) (3 si 

t=l 


Equations (3.5) and (3.4) can be put into a state space representation by defining 
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Modal 


Mode shape 

values at 


Mode # 

frequency 
rad /sec 

x=0.5 ft 

x=2.5 ft 

x=6.0 ft 

x=9.5 ft 

1 

0 

1.000 

1.000 

1.000 

1.000 

2 


-1.590 

-1.010 

0.000 

1.010 

3 

11.418 

-1.600 

-0.123 

1.210 

-0.123 

4 



-0.876 

0.000 

0.876 

5 

61.258 


-1.300 

1.410 

-1.300 

6 

100.900 


1.090 

0.000 

-1.090 

7 

150.185 


0.356 

1.400 

0.356 

8 

209.004 


-0.553 

0.000 

0.553 


Table 3-1: Modal Frequencies and Normalized 

Mode Shapes (taken from reference [4]) 


the control vector 


u T = [ Ui u 2 u 3 u 4 ] 


the measurement vector 


y T = [ yi ?3 y* I 

and the state vector 


x T = [ip l ip 2 $ 2 . . . ] 

Then, the state space representation of the beam model is 


i(l) = Ax(l) + Bu(l) (3.6) 

y(>) = c x(t) 


where 
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0 

« 2 (f 1 ) o ■ 

<y e i) 

0 

■M' 2 ) 

0 

<M< 2 ) o ■ 

• • <M f 2> 

0 


0 

Hh) 0 ■ 

■ • <W f 3) 

0 

) 

0 

*,(<„) 0 • 


0 



where — 0.5 ft., c 2 = 2.5 ft., * 3 = 6.0 ft. and e 4 = 9.5 ft., and N is the number 
of modes included in the model. 

The block diagonal structure of matrix A results from the zero coupling 
between the modes, which in turn, is due to the lack of damping in the beam. The 
structure of matrices B and C correspond to having force actuators and deflection 
sensors respectively. It is important to note that, any flexible structure with 
negligible damping having force actuators and displacement sensors, will be 
described by a state model with matrices A, B and C of the same structure as the 
ones corresponding to the beam. Therefore, most of the results obtained in this 
thesis can be applied to most flexible space structures. 

3.2 Filter Design for Actuator Failure Events 

The actuator failure model v/as described in section 2.1 and is given by 
equation (2.5). The event vector f • associated with the failure of actuator j is equal 
to by , the jth column of matrix B: 

fy=by ;=1, 2, 3, 4 


Let’s first consider the design of a failure detection filter to detect and 



identify the four actuator failure events. Applying the definition of output 
separability given in section 2.3.2, it was found that u^ = 1 and that the events { 
fj, f 2 , f 3 , f 4 } are output separable if and only if the number of modeled modes N is 
equal to or larger than four: N > 4. The detection generator corresponding to 
these events was found to be equal to f • , and the detection order equal to two for 
any N > 1: 


g. = f . = b • 

; j 


„.= 2 


3,4 


The number of completion spaces is (rank C - r) . Then, since rank C = 4 
and r == 4, there are no completion spaces, and, therefore the group detection order 
is equal to 2 X N, the number of filter states. Moreover, since the sum of the 
individual detection orders v- is 8, it can be concluded from the group detection 
theorem that the events { fj, f,,, f 3 , f 4 } are mutually detectable if N = 4. If N > 
4, the set of four actuator failure events win not be mutually detectable, but the 
set of any combination of up to three of these events will be. 

A failure detection filter was designed to detect the failures of actuators 
number 3 and 4. The design was made using the algorithms of Chapter 2 and was 
based on a model of the beam containing the first 7 modes, resulting in a 14 state 
filter. This design produced a state space partitioned in this way: 


R 14 = Ea © R, © e, © c 2 


where and are the two dimensional detection spaces corresponding to 
actuators number 3 and 4 respectively, and C t and C 2 are six and four dimensional 
completion spaces. In this first design, the 14 filter poles were rather arbitrarily 
placed at -10 rad/sec. This filter is named Detection Filter #1. 
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Detection Filter #1: 

# of modes = 7 

failure events = Actuators 3 and 4 
Poles: 

4 detection space poles at -10 rad/sec 
10 completion space poles at -10 rad/sec 

3.3 Computer Simulation 

In all the computer simulations, the system and filter start with zero initial 

conditions at time 0 sec.. The controls applied by the actuators are randomly- 

chosen command forces uniformly distributed between -1 and 1 pound, and held 

constant during 1/32 seconds. The actuator failure occurs at t = 1 sec. and the 

failure mode chosen is complete failure, that is \ij(t) = 0 or equivalently, rj(l) = 

u c (/) in equation (2.4). 

; 

The dynamics of the system and filter are simulated using a fifth order 
Runge-Kutta integration routine. The step size is chosen automatically by this 
routine, and is always smaller than 1/192 seconds. Therefore, since the period of 
the highest frequency mode being simulated is 1/33.2 seconds, the continuity of the 
system and filter dynamics is preserved. 

The detection filter #1, described in the previous section, is first tested with 
no model mismatch (the evaluation model of the beam consists of the first 7 
structural modes.) When the actuators are operating nominally, the output error is 
zero. When actuator 3 fails, the output error lies in the direction CA U 3 f 3 = 
CAb 3 , and when actuator 4 fails, it lies in the direction CA U 4 f 4 = CAb^ (since u 3 


r 
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= u 4 = 1.) One way of reading this directionality information from the output, in 
order to identify the failed actuator, is by transforming the output space: 

E(t) = Re{t) (3.7) 

where e(t) is the four dimensional output error vector, 
e(') = [e 3 E 4 ] T 

is the two dimensional transformed output error vector, and 
R = [(CF) t CF]' 1 (CF) t 
is the pseudoinverse (F is defined in section 2.3.2). 

With no model mismatch and before the failure occurs, the output error f(/) 
and transformed output error E(/) are both zero. When actuator 3 fails e(/) = 
CAb 3 K(t) , where k ( t) is a scalar time function. Therefore, 

E(f) =Rf(l) = RCAb 3 k( /) 

and since CAb 3 = CF 
then 

’ 1 ' 

E(<) = [(CF) t CF } A (CF) t CF K (t) 

L° 



and 


F 
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E(<) = 



(3.8) 


Hence, when actuator 3 fails, E 3 becomes nonzero and E 4 remains zero. In a 
similar way, it can be shown, when actuator 4 fails, E 4 becomes nonzero and E 3 
remains zero. 



Figure 3-2: Detection Filter #1. No model mismatch. 

Actuator 3 failure at T = 1 sec. 

These results were verified with simulations. Figure 3-2 shows the 
transformed output error (E 3 and E 4 ) resulting from the simulation of the beam 
and detection filter #1 with no model mismatch, and with the failure of actuator 
number 3 at T = 1 sec.. The transformed output error is zero until T = 1 sec., at 
which point E 3 starts growing, clearly indicating the failure of actuator number 3. 
Figure 3-3 shows the results of the simulation with no model mismatch and with 


failure of actuator number 4 at T = 1 sec.. 


a 



Figure 3-3: Detection Filter #1. No model mismatch. 

Actuator 4 failure at T = 1 sec. 

With no model mismatch, the detection filter #1 performed as expected, 
producing evident signatures at the output indicating the failure of actuator 
number 3 or 4. In the next chapter, the performance of this filter in the presence of 
model misrr ^tch will be analyzed. 



F 


-46- 


Chapter 4 

Failure Detection with Model Mismatch 

The best way of introducing the problem of failure detection in the presence 
of model mismatch, is by showing the results of the simulations of Detection Filter 
#1, introduced in the previous chapter, with mismatched evaluation models of the 
beam. Model mismatch occurs in two forms: unmodeled dynamics resulting from 
the truncation of the system model, and parameter errors in that portion of the 
system dynamics that is being modeled. 

In the flexible beam, unmodeled dynamics take the form of unmodeled high 
frequency structural modes. Since the model of the flexible beam is essentially 
infinitely dimensioned, the design of the detection filter is based on a reduced order 
model of the beam. Then, those high frequency structural modes that have been 
truncated from the model are the unmodeled modes. 

To test the performance of the detection filter in the presence of unmodeled 
modes, simulations are performed in which the evaluation model of the beam is of 
higher o^der than the model used in the filter design. For example, Detection Filter 
#1, based on a 7 mode model, is tested with an 8 mode evaluation model cf the 
beam. Hence, in this simulation the eighth structural mode, with a frequency of 
209.004 rad/sec., is the unmodeled mode. All the system modes, including the 
unmodeled one, are given zero initial conditions at T = 0 sec., and since the real 
physical beam has some damping, this is not an unrealistic assumption. Also, in 
order to highlight the effect of model mismatch, neither sensor noise nor 
disturbance is simulated throughout this thesis. 
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The result of the simulation of Detection Filter #1, with an 8 mode 
evaluation model of the system, and with failure of actuator 4 at T = 1 sec., is 
shown in Figure 4-1. The residual due to the unmodeled dynamics (the eighth 
structural mode) completely obscures the signature of the actuator failure (see also 
Fig. 3-3). Strikingly, this residual has a low frequency character rather than the 
frequency of the unmodeled bending mode. 


CJ 



Figure 4-1: Detection Filter #1. Eighth Mode is Unmodeled. 

Actuator 4 Failure at T = 1 sec. 

The parameter errors in the model of the flexible beam are the errors in the 
calculated values of the frequencies and shapes of the structural modes. That is, 
they are the errors in the figures of Table 3-1. Since the model frequencies and 
shapes depend on the physical characeristics of the beam, the errors in these 
quantities are correlated. However, when testing the performance of detection 
filters in the presence of parameter uncertainties percentage errors are introduced 
in the parameters of the evaluation model in a random manner with a uniform 
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distribution. For example, when considering parameter uncertainties of 5%, the 
percentage errors in the natural frequencies and shapes of the modes are uniformly 
distributed between -5% and +5%. 

The errors in each of the matrices A, B and C, resulting from the errors in 
the knowledge of the modal frequencies and shapes, affect the performance of the 
detection filter in different ways. Therefore, in order to gain a better 
understanding of the problem, the detection filters are tested with parameter 
uncertainties in only one of the matrices A, B or C. In all cases, the detection filters 
are designed using the values of Table 3-1, and tested with parameter uncertainties 
by introducing errors in the matrices A, B and/or C of the evaluation model of the 
beam. 

The result of the simulation of F etection Filter #1 with 0.05 % parameter 
uncertainties in the system matrix A (modal frequencies), and with failure of 
actuator 4 at T = 1 sec., is shown in Fig. 4-2. Notice that the noise due to the 
parameter errors is larger than the failure signature, even though the level of 
parameter uncertainties is extremely low. 

Figure 4-3 shows that the result of the simulation of Detection Filter #1 with 
0.05% parameter uncertainties in the control matrix B, and with failure of actuator 
4 at T = 1 sec.. The result of this simulation is even more surprising than the 
preceding ones. The residual due to errors smaller than 0.05% in the control 
matrix B, is of the same magnitude as the residual due to the complete failure of 
actuator 4 ! 

Finally, Fig. 4-4 shows the result of the simulation of Detection Filter #1 
with 0.05% parameter uncertainties in the measurement matrix C, and with failure 
of actuator 4 at T = 1 sec.. It can be seen that the residual due to the parameter 
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Figure 4-2: Detection Filter #1. Parameter Uncertainties in Matrix A = 0.05 

Actuator 4 Failure at T = 1 sec. 


PI 



Figure 4-3: Detection Filter #1. Parameter Uncertainties in Matrix B = 0.05*^. 

Actuator 4 Failure at T = 1 sec. 

errors in matrix C, is two orders of magnitude larger than the failure signatures. In 
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addition, notice the very low frequency character of this residual. It seems so far 
that the errors in the measurement matrix C are the most critical ones. 



Figure 4-4: Detection Filter #1. Parameter Uncertainties in Matrix C = 0.05 r c. 

Actuator 4 Failure at T = 1 sec. 

Figures 4-1 through 4-4 dramatically illustrate the effects of model mismatch 
on the performance of the detection filter. Even with extremely low model 
uncertainties and high frequency unmodeled modes, the residual due to model 
errors makes the failure signature indistinguishable. 

The objective in this chapter is to improve the visibility of the failure 
signature in the presence of mode! errors. To meet this objective, first the different 
dynamics associated with this problem are described. A frequency domain analysis 
that will lead to a more intelligent choice of the filter pole locations then follows. 
Finally, an approach to the design of robust detection filters based on low order 
models of the beam and output filtering, is presented. 


-51- 


4.1 Detection Filter Dynamics 

In order to understand the reason for the poor results shown in Fig. 4-1 
through 4-4, different aspects of the dynamics of detection filters are discussed in 
this section. Most of these results, however, are only valid when the plant is an 
undamped flexible structure, with a state model representation given by the 
matrices A, B and C described in Chapter 3. 

4.1.1 Structural Mode Dynamics 

The modal amplitude rp^t) of the ith mode, is given by equation (4.4) where 
the term to the right of the equal sign will be referred to as z -. The block diagram 
corresponding to this equation is shown in Fig. 4-5 



Figure 4-5: Structural Mode Dynamics 


The transfer function from z i to rp- is 
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(4.1) 


and the Bode plot of the frequency response T^;w) is given in Fig. 4-6. 



Figure 4-6: Modal Frequency Response 

Note in Fig. 4-6, that the frequency response of the mode approaches l/~'" as 
the frequency tends to zero. 


4.1.2 Detection Space Dynamics 

Given a detection filter configured to detect the failure event f, Beard showed 
that when event f occurs, the Laplace transform of the output error is: 

£ ( s ) = CA“ f H(s) r)( s) (4.2) 

v-i ■ • - +o - 

s'f + p sV‘ + . . . + P, 


H(s) = 
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where the zeros of the denominator polynomial are the poles of the detection filter 
associated with the detection space R^ , and the o^’s were defined in equation 
(2.19) . 

From equation (4.2), it can be concluded that the transfer from rj( s) to c(s) 
depends only on the dynamics associated with the detection space R f . This result 
is not surprising, since R f is an invariant space and f c R f . Also note that while the 
designer has complete control over the denominator of H(s), by choosing the 
detection space poles, he cannot alter the numerator. 

These results can be specialized to the problem of detecting the failure of an 
actuator in the flexible beam. It was found, for actuator failure events, that u = 1. 
i/j = 2 and g — f (equivalently Oj = 1). In particular, for Detection Filter #1 
(with all detection space poles at -10 rad/sec.) and complete failure of actuator 4 
(tj(<) = u 4 (/) ) equation (4.2) becomes 

f(s) = CAb 4 H 4 (s) u 4 (s) (4.3) 

where 



1 

s 2 + 20 s -I- 100 


and 


Moreover, since the transformed output error is 
E(s) = R t(s) = R CAb 4 H 4 (s) u 4 (s) 
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R CAb 4 = 


0 

1 


then 

E 4 (s) — H 4 (s) u 4 (s) (4.4) 

The Bode plot of the frequency response H 4 (yo;), corresponding to Detection 
Filter #1, is given in Fig. 4-7. 



Figure 4-7: Detection Space Transfer 


4.1.3 Model Error Dynamics 

In this subsection, models for the different types of model mismatch are 
developed. These models are the basis for the frequency response analysis 
performed in the next section, that will help explain the very large residuals 
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produced by ” v ery small” modeling errors. 

In all the following cases of model mismatch, the state and output equations 
of the filter are the same as before, and for convenience are repeated here: 

i(l) = A x(() + B u c (/) + D [ y(l) - y( <) ] (4.5) 

y(f) = Ci(() 


where the matrices A, B and C were defined in the previous chapter for N modes. 
First, the case in which the system has more than N modes (this will always be the 
case for flexible structures) is analyzed. Then, three cases are analyzed in which 
only one of the matrices A, B or C contain errors and the other two are exact. 
Finally, the real situation in which the three matrices contain errors is considered. 

4.1.3. 1 Unmodeled Modes 

Suppose that N is the number of structural modes that are modeled in the 
design of the detection filter. Then, the state and output equations for the beam 
can be written as 



u 


where the matrices A, B and C correspond to the modeled part of the system 


dynamics, and A u , B u and C u correspond to the unmodeled part. The state and 
output equations for the filter are exactly the same as before (Equation (4.5)). 
Therefore, the new state and output error equations are: 

t(«) = (A-DC)«(0-DC.*.(0 

<(()=C.(l) + C,x,(l) (4.6) 

where 

*,(') = A„x„(f) + B„u(<) 

The poles of the filter are still given by the eigenvalues of (A - DC), and therefore, 
the presence of unmodeled modes does not affect the stability of the filter. They 
do, however, impair the detecton capabilities of the filter, by generating a residual 
that obscures the failure signature. 

For the sake of clarity, equation (4.6) is written as 

CO 

4(4) = (A- DC) e(()- D £ (4.7) 

»=N+1 


where 
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and where ip- (f) is the amplitude of an unmodeled mode. 



Figure 4-8: Unmodeled Modes Error Dynamics 

Figure 4-8 shows a block diagram representation of equation (4.7). The actuators 
drive the unmodeled modes in the same way shown in the block diagram of Fig. 
4*5. The unmodeled modes amplitudes ^n+I’ ^N+°’ ' ' ' ’^en generate an output 
residual through two different paths: a direct path through the vectors c^, and an 
indirect path by driving the state error equations through the "control vector" Dc, 
Note that while the designer has no control over the direct path , he has "some” 
control on the indirect path by adjusting the filter pole locations with the gain 
matrix D. However, since D must be a detector gain, and is in both the feedback 
loop (A - DC) and the control vector Dc { , this is not a trivial problem. 

4. 1.3.2 Parameter Errors in A 

The parameters of the system matrix A are the modal natural frequencies. 
The true modal frequencies are called oj-, and the true system matrix is called A*. 







r 
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Hence, if u i and A are the modal frequencies and the system matrix used in the 
filter design, define: 

o *2 n 

A ojrz=u i - ujf (4.8) 

AA = A* - A 

and, therefore, 



Hence, the state and output equations for the system are now: 
i(() = (A + AA) x(f) + B u(<) 

y(f) = C x(/) (4.9) 


Combining these equations with the equations of the filter (Eq. (4.5)), the new state 
and ouiput error equations are: 
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i{t) = { A -DC) e(<) + AA x(f) 

t(/) = Ce(() (4.10) 


These equations can also be written as: 


N 

«(l) = (A-DC)e(f)+£Aa,.VK(<) 
((f) = C e(t) 


(4,11) 


where 


0 


Aa i= Au i ®N(2 X i) 


0 

Au/“ 


0 


0 


The block diagram representing the equation (4.11) is illustrated in Fig. 4-9. 

From Eq. (4.11) or Fig. 4-9, notice that the modal amplitudes 0. drive the 
state error through the "control vector” Aa^ , thus generating a residual at the 
output. This error model is similar to the error model for actuator failures, but 
with 0. and Aa. in place of u ■ and b^. 
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Figure 4-9: Error Dynamics for Model Error in A 
4. 1.3. 3 Parameter Errors in B 

The parameters of the control matrix B are the values of the mode shapes at 
the actuator locations. The true control matrix is B , and the control matrix used 
in the filter design is B. Define the difference between these two matrices as: 

AB = B* - B 

where 

AB = [ Abj Ab 2 Ab 3 Ab 4 ] 

Therefore, the system state and output equations are now: 


*(/) = Ax(l) + (B + AB) u {t) 


(4.12) 







y(0 = c *(/) 
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and, combining them with the filter equations, get the new state and output error 
equations: 


«(() = (A- DC) «(/) + £>b,. 

<(<) = C«(0 


(4.13) 


These equations are represented in a block diagram form in Fig. 4-10. 



Figure 4-10: Error Dynamics for Model Error in B. 

From Eq. (4.13) and Fig. 4-10, it can be noted that this error model is the 
same as the error model for actuator failures but with Ab t - instead of b ( . 
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4. 1.3.4 Parameter Errors in C 

The parameters of the measurement matrix C are the values of the mode 
shapes at the sensor locations. Define C* as the true measurement matrix. The 
difference between C* and C (the measurement matrix used in the filter design), is 
defined as: 

AC = C* - C 

where 

AC = [ ACj AC 2 . . . AC n ] 

Then, the state and output system equations are: 

i(t) = Ax(f) + Bu(i) (4.14) 

y(f) = ( C + AC ) x(i) 

and the new state and output error equations are: 

N 

4(0 = (A - DC) e(f)-D£ Ac. V-, (i) (4.15) 

4=1 

N 

£ W = Ce(l) + ^Ac 
4=1 

In Fig. 4-11, the block diagram representation of Eq. 4-15 is shown. It is important 
to note that this error model is similar to the error model for unmodeled modes 
(Fig. 4-8). The main difference is that while the state error in Eq. (4.15) is driven 
by the amplitude of the modeled structural modes, the state error in Eq. (4.7) is 
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driven by the amplitude of the unmodeled ones. The other difference is in the 
vector multiplying the modal amplitudes: in Eq. (4.15) this vector is AC^, and in 
Eq. (4.7) is 



Figure 4-11: Error Dynamics for Model Error in C 

4. 1.3.5 Parameter Errors in A, B and C 

In the real problem, all the matrices A, B and C contain errors. This problem 
can be analyzed by dividing it into the three cases explained in the three previous 
subsections, and then adding the output residuals obtained in each case. However, 
in order for the superposition principle to apply, a modification must be introduced 
when dividing the problem into the three cases illustrated in Figures 4 9, 4-10 and 

4-11: the modal amplitudes tp i in Fig. 4-9 and 4-11 must be generated by the true 

* * 

beam model, with the true system matrix A and the true control matrix B 

instead of the matrices A and B used in the filter design. 

The modal amplitudes obtained with A and B*, however, are not different in 
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character than the ones obtained with A and B (for practical values of parameter 
errors). Moreover, the reasons for the extremely large residuals when there are 
parameter errors in the model, lie on the transfers from the modal amplitudes i> i to 
the output t(t) in Fig. 4-9 and Fig. 4-11, and on the transfers from the control 
inputs u . to the output f(f) in Fig. 4-10. Therefore, the problem of parameter 
errors in the matrices A, B and C can be analyzed by dividing the problem into 
three cases in which only one of the matrices contains errors and the other two are 
exact. This is the approach taken in this thesis. 

4.2 Frequency Response Analysis 

The reasons for the poor performance exhibited by Detecton Filter #1 in the 
presence of model mismatch, can be explained by conducting a frequency response 
analysis, based on the error models developed in the previous section. 

From Eq. (4.1), it can be concluded that the higher the natural frequency of 
the structural mode, the lower the amplitude transfer at low frequencies. 
Therefore, since the unmodeled mode (the eight structural mode) is a very high 
frequency mode, its amplitude has very small low frequency components, and most 
of its energy is at its natural frequency. However, Fig. 4-1 shows, very 
surprisingly, that the unmodeled mode produced a residual with a low frequency 
character rather than the high frequency of the unmodeled bending mode. 

The reason for this residual behavior can be explained in terms of the 
frequency response functions for the transfer from the amplitude of the eight mode 
V’g to the residuals E 3 and E 4 (see Fig. 4-8). These frequency response functions 
are shown in Fig. 4-12, where the extremely high amplitude transfer at low 
frequencies can be noted. This very large gain at low frequencies amplifies the 
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small low frequency components of the unmodeled mode excitation, giving the 
residuals of Fig. 4-1 a low frequency character. 

A simple calculation can show that the low frequency gains in Fig. 4-i2 result 
in an unfavorable signal to noise ratio. This calculation is only an approximation 
and is based on the assumption that actuator 4, the one that fails at T — 1 sec., is 
the only active actuator. Under this assumption, the low frequency gain for the 
transfer from u 4 to ^ g can be computed as: 

si b « 1102 , 

— - = ~ = = 2.522 X 10' 5 for u) « 0 

|u 4 | 43683 


where b g = <p & (e 4 ) / M (see Fig. 4-5). The low frequency gain for the transfer from 
V> g to E 4 can be obtained from Fig. 4-12: 



for u) ^ 0 


Therefore, the ratio |EJ / juj is calculated as: 


|e 4 ! 


H = — = 2.522 X 10* 5 X 2070 = 0.052 for u % 0 

" l« 4 l 


This ratio corresponds to the noise present at the residual E 4 , due to the 
unmodeled mode. The ratio corresponding to the signature produced by the failure 
of actuator 4 is, from equations (4.3) and (4.4): 
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H 


s 



for u 0 


Therefore, the low frequency signal to noise ratio at E 4 is H s /H n = 1/5.2. This 
unfavorable signal to noise ratio explains why low frequency noise obscured 
completely the failure signature in Fig. 4-1. 

In the real experiment, the four actuators are active. In this case, it is 
reasonable to expect the modal amplitude V’g to be larger than the one resulting 
from only one actuator. Consequently, a larger noise at is also expected. The 
failure signature, however, remains the same because it depends only on actuator I. 
Therefore, it is likely in the real experiment, for the signal to noise ratio to be 
worse than 1/5.2. 

The results shown in Fig. 4-12 were very unexpected. A detection filter is a 
linear observer whose gain matrix D has the additional constraint, besides making 
the filter stable, of restricting the output error due to a particular failure to a 
single line in the output space. The typical behavior of a linear observer is: 1) to 
follow the measurement noise when its frequency lies inside the bandwidth of the 
filter, and 2) not to respond to the measurement noise when its frequency is higher 
than the filter bandwidth, thus resulting in large residuals. Translated into the 
frequency domain, this means that the frequency response of the transfer from the 
measurement noise to the residual is expected to be very small within the 
bandwidth of the filter, and then grow to a certain value at higher frequencies. 
Detection Filter #1 behaved as expected at frequencies higher than 10 rad/sec by 
responding moderately to the unmodeled mode noise. However, it behaved quite 
differently than expected at low frequencies, by actually amplifying the unmodeled 
mode noise, as can be sec-n in Fig. 4-12. It seems that the directionality constraint 
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Figure 4-12: Detection Filter #1. Unmodeled Mode Transfer 
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con.bined with the dynamics of the flexible beam, resulted in a detection filter that 
is a very poor observer when all of its poles are placed at -10 rad/sec. 

The results of testing Detection Filter #1 with 0.05% parameter 
uncertainties in system matrix A, can be seen in Fig. 4-2. Note that the residual 
that obscures the failure signature has the frequency of the first bending mode. 
The error dynamics for the case of parameter errors in A is shown in Fig. 4-9. The 
vectors and Aa 0 are zero because they correspond to the rigid modes, whose 
frequencies are known to be exactly zero. Of the bending modes, the first one, with 
a frequency of 11 rad/sec, comes the closest to lying within the bandwidth of the 
filter. Therefore, since the contributions from the higher bending modes are 
filtered out, the residuals have the frequency of the first bending mode. Figure 
4-13 shows the frequency response of the transfer from V> 3 to the residuals E 3 and 

£ 4 ‘ 

When testing Detection Filter #1 with parameter uncertainties, the most 
surprising result was obtained when considering parameter errors in matrix B. The 
residuals produced by parameter errors smaller than 0 05% in the control matrix B, 
were of the same magnitude as the signature resulting from the complete failure of 
actuator 4 (See Fig. 4-3). 

As can be seen in Fig. 4-10, the reason for these large residuals must lie in the 
transfers from the control imputs u $ to the residuals E 3 and E^. The frequency 
response functions for the transfer from the control u 4 to the residuals E 3 and E 4 
are shown in Fig. 4-14. The values of these functions at low frequency are 0.0065 
for E 3 and 0.011 for E 4 . Therefore, since the low frequency gain corresponding to 
failure of actuator 4 is 0.01, the residuals and the failure signature are of the same 
magnitude. 
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Figure 4-13: Detection Filter #1. Parameter Errors in A = 0.05%. 
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It is interesting to note that even though the elements of the vector Ab 4 are 
smaller than 0.05% of the respective elements in b 4 , the amplitude transfer from u 4 
to E 4 is actually a little larger with Ab 4 than with b 4 as the control vector in F le;. 
4-10. However, this would not be the case if Ab 4 lay in the same direction in the 
state space as b 4 . For example, if the error in each element of b 4 were +0.05%, 
then Ab 4 = 0,0005 X b 4 and the gain corresponding to Ab 4 would be 0.05% of 
the gain corresponding to b 4 (i.e., 0.0005 X 0.01 = 5 X 10' 6 ). 

Since the dynamics of both detection spaces are the same, then the large low 
frequency gains that result when Ab 4 does not lie in the same direction as b 4 must 
be due to the dynamics of the completion spaces. 

The test of Detection Filter #1 with 0.05% parameter uncertainties in the 
measurement matrix C, produced the results shown in Fig. 4-4. The residuals have 
a very low frequency character and appear to be unbounded. 

The error model for the case of parameter errors in C, is shown ir« Fig. 4-11. 
As mentioned earlier, it is similar to the error model for unmodeled modes. Figure 
4-15 shows the frequency response functions for the transfer from the modal 
amplitude rp^ to the residuals Eg and E 4 - The low frequency gains are almost three 
orders of magnitude smaller than those corresponding to the transfer from the 
unmodeled mode rp^ to the residuals E 3 and E 4 . The reason for this is that the 
vector Acj is much smaller than the vector c 8 - 

The beam is being driven by the actuators in an open loop manner. 
Therefore, the amplitude of the rigid modes tp^ and ip n} since their poles are at the 
origin, start to increase and become unbounded. That is, the beam moves away 
from the reference and starts to rotate. Therefore, even though the low frequency 
transfer from ip l and ip 2 to E is relatively small, since the modal amplitudes xp^ and 
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Figure 4-14: Detection Filter #1. Parameter Errors in B 

Transfer from u 4 to E, and E 4 . 








if; 0 are growing with no bounds, eventually and V> 2 become large enough to 
generate residuals that obscure the failure signature. This explains the very low 
frequency and the unbounded nature of the residuals in Fig. 4-4. 

To summarize, the reasons for the poor performance of Detection Filter #1 in 
the presence of model mismatch were explained in terms of frequency response 
functions (Fig. 4-12 through 4-15). It was found that these frequency responses had 
features that were very unexpected, like very large low frequency gains. 

It is believed that the "strange” transfer properties of Detection Filter #1 are 
the result of applying the directionality constraints to the dynamics of flexible 
structures. The question now is whether a different choice of filter pole locations 
would improve these transfer properties and the performance of the filter. 

It is reasonable to expect that by increasing the bandwidth of the filter, by 
moving the filter poles further to the left in the s - plane, the filter's ability to 
track the measurement vector y(t) will improve, thus resulting in smaller error 
residuals. This would be true for both detection space and completion space poles. 
However, since moving the detection space poles to the left also reduces the size of 
the failure signature (Eq. (4.2)), it is not clea*- whether this would result in a better 
signal to noise ratio. Conversely, since the failure signature does not depend on the 
completion space poles, moving these poles to the left could result in a better signal 
to noise ratio. 

To test this hypothesis, a detection filter was designed with the same 
detection space poles, but with the completion space poles moved to -100 rad/sec. 

Detection Filter #2: 

# of modes = 7 


failure events = Actuators 3 and 4 
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Figure 4-15: Detection Filter #1. Parameter Errors in C = O.OS^. 

Transfer from to E 3 and E^. 
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Poles: 

4 detection space poles at -10 rad/sec 
10 completion space poles at -100 rad/sec 

Since the location of the detection space poles was not changed, the 
signatures for the failure of actuators 3 and 4 remained the same as before, and are 
shown in Fig 3-2 and 3-3. 

For the remainder of this section, Detection Filter #2 will be subject to the 
same frequency response analysis and will be simulated under the same conditions 
as Detection Filter #1. 

Figure 4-16 shows the frequency response functions for the transfer from the 
unmodeled mode to the residuals E 3 and E 4 , for Detection Filter #2. Note the 
dramatic decrease in the low frequency gain. . 

The improved transfer Droperties of Detection Filter #2 suggest that it 
should perform much better than Filter #1 in the presence of unmodeled dynamics. 
This is shown to be the case in Fig. 4-17, where Detection Filter #2 is tested in the 
presence of the unmodeled eight mode. The residuals no longer have a low 
frequency character, and the frequency of the bending mode is evident. More 
importantly, the signature of the failure of actuator 4 at 1 sec. is clearly visible. 

Figure 4-18 shows that there has also been a dramatic improvement in the 
transfer corresponding to the case of parameter errors in 

The results of the simulation of Detection Filter #2 with 0.05% parameter 
uncertainties in A, are shown in Fig. 4-19. The large improvements in the transfer 
properties are reflected in these results, where the failure signature now is clearly 
visible. The bias in the residuals has vanished and the frequency of the first 
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Figure 4-16: Detection Filter #2. Unmodeled Mode Transfer 



-76- 



Figure 4-17: Detection Filter #2. Eighth Mode is Unmodeled. 

Actuator 4 Failure at T = 1 sec. 

bending mode is barely noticeable. Now, the frequency of the second bending mode 
becomes evident. 

Notice that the errors in the modal frequencies of the third, fourth and fifth 
bending modes do not produce any noticeable effects in the residuals. The reason 
for this desirable result is that the frequencies of these bending modes lie outside 
the bandwidth of the filter, and therefore, they are filtered out. Since the 
uncertainty in the knowledge of the modal frequencies increases the higher these 
frequencies are, this result is very important. 

Figure 4-20 shows the residuals produced by the simulation of Detection 
Filter #2 with \% parameter uncertainty in matrix A. The failure signature is 
completely obscured by the residuals, and therefore, it can be concluded that the 
improved filter cannot perform adequately with uncertainties equal or larger than 
\% in the knowledge of the modal frequencies. 
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Figure 4-19: Detection Filter #2. Parameter Errors in Matrix A = 0 OS' 7 ©. 

Actuator 4 Failure at T = 1 sec. 



Figure 4-20: Detection Filter #2. Parameter Errors in Matrix A = l r o. 

Actuator 4 Failure at T = 1 sec. 
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Detection Filter #2 also showed a remarkable improvement in performance 
when tested in the presence of parameter uncertainties in B. The frequency 
response functions related to this case are shown in Fig. 4-21. The low frequency 
gains are two orders of magnitude smaller than the gain corresponding to Detection 
Filter # 1 . 

The results of the simulation of Detection Filter #2 with 0.05% parameter 
uncertainties in B, can be seen in Fig. 4-22. The improvement is remarkable. The 
residuals due to model mismatch are barely noticeable. However, the improved 
filter is still not good enough to detect failures in the presence of more realistic 
levels of parameter uncertainties. The filter was tested with 5% parameter 
uncertainties in B, and the signature was completely obscured by the residual (Fig. 
4-23). 

Finally, Detection Filter #2 was tested with parameter uncertainties in the 
measurement matrix C. Its transfer properties are shown in Fig. 4-24, where a large 
improvement can also be noted. The simulation, shown in Fig. 4-25 gave much 
better results as well. However, the improvements were not large enough to make 
the failure signature clearly visible. 

To summarize, dramatic improvements have been accomplished by moving 
the completion space poles further to the left in the s - plane. The largest 
improvements were obtained in the performance of the filter in the presence of 
unmodeled modes and parameter errors in B. The improved filter is able to detect 
failures with almost 5 % parameter uncertainties in B. 

Even though the improvements in the performance of the filter with model 
errors in A and C were also remarkable, the filter would not tolerate parameter 
errors in A of up to 1%, and in the case of system matrix C, the improvements 
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Figure 4-23: Detection Filter #2. Parameter Errors in Matrix B = 5%. 

Actuator 4 Failure at T — 1 sec. 
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Figure 4-25: Detection Filter #2. Parameter Errors in C = O.OS^. 

Actuator 4 Failure at T = 1 sec. 

were not enough to allow successful detection with 0.05% parameter uncertainties. 

However, the reasons for this disappointing performance may lie in part on 
the somewhat unrealistic conditions under which the filter was tested. In the case 
of parameter errors in A, the state error was driven mainly by the very large 
components of the modal amplitudes at their natural frequencies resulting from the 
lack of damping in the beam. The real beam, however, has some damping and 
therefore these components would be considerably smaller. In addition, all the 
modal amplitudes are kept small by the action of the control system during its 
nominal operation. In the case of parameter errors in C, the state error was driven 
by the unbounded amplitudes of the rigid modes corresponding to the drifting and 
rotation of the beam. In the real problem, however, the control system would 
prevent the beam from drifting and rotating. 


It is also worth noting that the performance of the detection filter with 


r 
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parameter errors in A could be improved by filtering the output. Figure 4-20 shows 
that the error residuals have a very large component at the frequency of the second 
bending mode (31.36 rad/sec.). Therefore, since the bandwidth of the failure 
signature is 10 rad/sec., the performance of Detection Filter #2 could be improved 
by filtering E 3 and E 4 with a bandwidth of 10 rad/sec. . 

Finally, something should be said about the placement of the detecton space 
poles. Unlike the case of the completion space poles, the choice of the locations of 
the detection space poles strongly depends on the frequency content of the control 
signals. It is clear that moving these poles to the left in the s - plane reduces the 
residuals due to model mismatch and reduces the failure signature as well. 
However, whether the signal to noise ratio improves or not will depend, among 
other things, on the frequency content of the control inputs. The same is true 
when considering the damping associated with these poles. Therefore, the choice of 
the location of the detection space poles should be made on the basis of more 
realistic simulations which include the control system along with the disturbances 
of the beam. Moreover, a consideration to be taken into account when choosing 
these poles is the response time of the detection filter. 

4.3 Reduced Order Detection i ilters 

In Chapter 3, it was shown that detection filters based on a model of the 
beam with more than 4 structural modes (the number of sensors) could not be 
configured to detect more than three actuator failure events. Moreover, the design 
of these filters was made very difficult by the presence of completion spaces, and 
the specification of the filter poles did not result in a unique gain matrix 
D. Another disadvantage of these filters, as shown by Detecton Filters #1 and #2, 
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was their very complex dynamics. Their robust properties varied dramatically with 
the location of the completion space poles, and even the improved filter (Detection 
Filter #2) could not tolerate realistic levels of parameter uncertainties. 

This section analyzes the design and performance of detection filters based on 
a model of the beam containing as many structural modes as the number of 
displacement sensors p, and with all of its poles placed in the same location. It will 
be shown that these filters do not have the disadvantages mentioned above. Also, 
ways of coping with the larger number of unmodeled modes will be discussed, along 
with the effects of using a reduced order model upon the detection and 
identification capabilities of the filter. 

4.3.1 Detection Filters with N = p 

Let’s consider the design of a detection filter based on a model of the beam 
with p structural modes, where p is the number of independent displacement 
measurements. That is, 

N = p = # of sensors — rank C 


Define the p X p matrix M as: 
M = [ Cj c 2 . . . c p ] 


where c ; is the column of matrix C corresponding to the amplitude of the ith 
structural mode. It will be assumed that rank M = p, and therefore, M is 
nonsingular. If this is not the case, the number of sensors used in the detection 
filter should be reduced until rank M = p. 
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In order to simplify both the filter design and the analysis that will follow, a 
transformation is performed on the measurement vector y: 

y + = M* 1 y 

then, 

y + = M* 1 C x 

and 

y + = C + x 

where 


C + = M' 1 C = 



0 0 

1 0 

0 0 

• * 

• • 

0 0 



From the structure of C + , it is easy to verify that 
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^1 

^2 


(4.16) 


Figure 4-26 shows the new block diagram of the failure detection filter which 
incorporates the measurement transformation M* 1 . Notice that the detection filter 
now uses the transformed measurement matrix C + . Also, it is easy to show that: 

D = D + M’ 1 
( = M c + 

The state and output error equations corresponding to the failure of the jth 
actuator are: 


«(() = ( A- D + C + ) e(i) + by i](t) 
< + (f) = C + e (i) 

where the vector h- has the structure: 



0 byj 0 h- 2 



(4.17) 


(4.18) 


Theorem 4.1 : 

The gain matrix D + that is a detector gain for all actuator failure events b ■ , 
and that assigns the filter poles as given by the roots of the equation 

0 = (s 2 £ u>p s -f- 
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2^ F 



Proof: 

First, calculate the matrix (A - D + C + ) as: 


A - D + C + = 


0 


0 


( 4 . 20 ) 


where the 2X2 matrix P is: 


P = 






o 


The poles of the filter are given by the roots of the characteristic equation of the 
matrix (A - D + C + ). That is, the roots of the equation 

0 = det(s I - A + D + C 4 ) 

Since the matrix (A - D + C + ) is block diagonal, using a result from linear algebra 
get: 
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det(s I - A + D + C + ) = ( det(s I - P) ) N = (s 2 + 2£w F s + u| ) N 
Therefore, the poles of the filter are given by the roots of the equation 
0 = (s 2 + 2£u;p s + Wp ) N 
as stated in the theorem. 

To prove that the gain matrix D + given in (4.19) is a detector gain for all 
actuator failure events, the transfer function vector T v (s) defined as 

€ + (s) = Tj(s) T](s) ( 4 . 21 ) 

Tj (s) = C + (s I - A + D + C + )* J b ; 

is calculated. It is easy to show that 


Q 


(s I - A + D + C + Y = 


Q 


o 


o 


Q 


where 


Q — ( s I * P}” 1 — 


s 2 -f 2£u; F s + 


-Wp s 2^a»p 

Then, by premultiplying by C + and postmultiplying by b^ obtain: 
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T ; (s) = 


s' + s + Wp 


b j2 


jN J 


where b.j, b y 2 , . . . , are the elements of the event vector by in (4.18). Finally, 
by defining the vector b+as: 




(4.22) 


equation (4.21) can be written 


€ + (s) = b + - 

J 


n(s) 


s' 4- 2^Wp s + Up 


(4.23) 


Then, it can be concluded that when the jth actuator fails, the output error c + 
remains in the single direction given by the vector b "t Therefore, since the event 
vector by represents any actuator failure, the matrix D + in (4.19) is a detector gain 
for all actuator failure events, and the proof is complete. 


V 
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As it was shown in Theorem 4.1, the design of the detection filters being considered 
in this section is trivial. The gain matrix D, in the filter implementation without 
the transformation, (Fig. 2-1), can be obtained as D = D + M' 1 . This gain matrix D 
is a detector gain for all actuator failure events as well, and the output error e 
remains in the direction CAby when the jth actuator fails. 

So far, two important advantages of detection filters with N = p over higher 
order filters have been shown: 1) these filters can be configured to detect all 

actuator failure events, and 2) the filter design is trivial. Now, their performance 
in the presence of model mismatch will be analyzed. 

A detection filter with N = p = 4 was designed using Theorem 4.1 with all 
of its poles placed at -10 rad/sec (i.e., <x>p — 10 rad/sec and £ = 1). This detection 
filter was named Detection Filter #3. 

Detection Filter #3: 

# of modes = # of sensors = 4 

failure events: all actuator failures 

Poles: 8 detection space poles at -10 rad/sec 

To be able to compare the performance of this filter with the performance of 
Detection Filters #1 and #2, the output c is transformed with the pseudo inverse 
matrix R to detect only failures of actuators 3 and 4, as was done with the previous 
filters. That is, the matrix R in Fig. 4-26 is the pseudo inverse of the matrix CF: 

R = [(CF) t CF l* 1 (CF ) T (4.24) 

where 


CF = [CAb, CAb 4 ] 
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Since the detection space poles of Detection Filter #3 are the same as in 
Filters #1 and #2, and since the transfer function from u 4 to E 4 (for complete 
failure of actuator 4) is still given by Eq.(4.4), the failure signature produced by the 
complete failure of actuator 4 is exactly the same as in Fig. 3-3. 

It is easy to see that Detection Filter #3 is going to perform much better in 
the presence of parameter errors in control matrix B than Filters #2 and #1. The 
vectors Ab ■ in Fig. 4-10 have the same structure as the actuator failure events 
in (4.18). Therefore, the vectors A by generate unidirectional error residuals 
governed by Eq. (4.3) with Aby and Uy in place of b 4 and u 4 . This means that the 
error residual generated by Aby , when compared with the signature produced by 
by , is of the same order as the percentage errors in by . In other words, there is no 
direction-dependent amplification of A by as in Filters #1 and #2. 

Figure 4-27 shows the frequency response functions for the transfer from u 4 to 
the residuals E 3 and E 4 , corresponding to 0.05?o parameter errors in B. The low 
frequency gains are two c. iers of magnitude smaller than the gains for Detection 
Filter #2 shown in Fig. 4-14. Therefore, Detection Filter #3 will perform much 
better in the presence of parameter errors in B than Filters #1 and #2. 

It is also reasonable to expect Detection Filter #3 to perform better in the 
presence of parameter errors in system matrix A, for the same reasons it performed 
better with parameter errors in B. The vector Aa . in Fig 4-9 also have the same 
structure as the actuator failure events by, and therefore also generate 
unidirectional output through the dynamics of the detection spaces. 

Figure 4-28 shows the response function of the transfer from the unmodeled 
eighth mode to the residuals E 3 and E 4 . The low frequency gains are one order of 
magnitude smaller than for Detection Filter #2 in Fig. 4-16. Therefore, Detection 
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Figure 4-27: Detection Filter #3. Parameter Errors in B = 0.05^. 
Transfer from to and E^. R = Pseudo Inverse. 
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Filter #3 will perform better in the presence of the unmodeled eighth mode than 
Filter #2. But, on the other hand, the fifth, sixth and seventh structural modes 
are unmodeled in Detection Filter #3, while they are modeled in Filter #2. In 
addition, these unmodeled modes have smaller natural frequencies than the eighth 
mode, and therefore, they are more excited. Hence, the problem of unmodeled 
modes is more serious in Detection Filter #3 than in Detection Filter #2. 

Part of the problem of unmodeled dynamics in Detection Filter #3 can be 
solved by output filtering. Since the natural frequencies of the unmodeled modes 
are higher than the bandwidth of the filter (i.e., 10 rad/sec), by filtering the output 
error ( with a bandwidth of 10 rad/sec, the large components of the error residuals, 
due to the natural frequencies of the unmodeled modes passing through the direct 
path in Fig. 4-8, can be eliminated without affecting the failure signature. 
Therefore, the elements of the output error vector E were filtered by second order 
filters with bandwidths of 10 rad/sec. . 

Figure 4-29 shows the results of the simulation of Detection Filter #3 with no 
model mismatch, with failure of actuator 4 at T = 1 sec. and with output filtering. 
Note that the failure signature is smailer than the failure signature produced with 
no output filtering in Fig. 3-3. This difference could have been avoided by choosing 
a slightly wider bandwidth for the output filters. 

Detection Filter #3 was then tested in the presence of unmodeled dynamics. 
Figure 4-30 shows the results of the simulation of Detection Filter #3 with an 8 
mode evaluation model of the beam, with failure of actuator 4 at T = 1 sec. and 
with output filtering. The enor residuals due to the four unmodeled modes (the 
fifth, sixth seventh and eighth structural modes) completely obscure the failure 
signature. From the low frequency character of these error residuals, it can be 
concluded that the output filter was successful in eliminating the high frequency 
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Figure 4-29: Detection Filter #3. No Model Mismatch. Actuator 4 Failure at 
T = 1 sec. . R = Pseudo Inverse. Output Filtering. 
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Figure 4-30: Detection Filter #3. Four Unmodeled Modes. Actuator 4 Failure at 
T = 1 sec. . R = Pseudo Inverse. Output Filtering. 
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components corresponding to tbe modal frequencies passing through the direct 
path. Then, the large error residuals are due to the low frequency components of 
the unmodeled mode amplitudes that lie within the bandwidth of the filter, and 
therefore, cannot be filtered out without affecting the failure signature. 

In the next two sections, ways of improving the performance of Detection 
Filter #3 in the presence of unmodeled modes will be introduced. 

4.3.2 Improved Output Transformation 

An impressive improvement in the performance of Detection Filter #3 in the 
presence of model mismatch can be obtained by simply using a different output 
transformation R. But before introducing the new transformation, an analysis of 
the dynamics of the detection filter described in Theorem 4.1 is in order. 

To simplify the analysis, only one unmodeled mode will be considered: the 
eighth structural mode. Then, the state and output error equations are 

i(i) = (A-D + C + ) e(l) - D + Cg *>,(() 

£+ (0 — C + e(f) + Cg V’g (*) (4.25) 

e(t) = M £ + (t) 

where 

C Z=M A c 8 

and Cg is the column of the measurement matrix C corresponding to the amplitude 
of the eighth structural mode. Due to the structure of matrices A, D + and C + , 
Eq.((4.25)) can be written as: 
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t.= 

4 = 

( A ,-d+c T )e r d + eg,.*-, 

c T •«+«»<*» 

i=l,2,. 

. .N 

(4.26) 

where 






e ,= 

Mr $i). 

c T = [l 0] 

0 1 ' 

-J] 0 

and where dtwas defined in (4.19), c g| - is the ith element of the vector and ftis 
the ith element of the output error vector c + . The two elements of the vector e i 
are the errors in the estimate of the amplitude of the ith structural mode and the 
error in the estimate of its rate of change. The block diagram corresponding to Eq. 
(4.26) is shown in Fig. 4 31. 

From Eq. (4.26) and Fig. 4-31, it can be concluded that the error dynamics of 
the detection filter described in Theorem 4.1 are decoupled. That is, the detection 
filter estimates the modal amplitude 0 and its rate of change ty for each structural 
mode separately. This decoupling is achieved in the following way. First, since the 
number of modes is equal to the number of independent displacement 
measurements, it is possible to calculate the modal amplitudes by transforming the 
measurement vector y to y 4 as can be seen in (4.16). Second, the structure of the 
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Figure 4-31: Decoupled Dynamics of Detection Filters with N = p. 

gain matrix D + in (4.19) applies the residual ip i - to the corresponding ith mode 
(see Figure 4-26). 

The decoupled dynamics of the detection filter described in Theorem 4.1 are 
very easy to analyze. As can be seen in Fig. 4-31, the unmodeled modes can be 
treated as disturbances in the decoupled loops corresponding to each of the 
modeled modes. Then, the transfer from V> 8 to c+can be easily obtained as: 

€t( s ) s 2 + J 

— = (4.27) 

c 8 i s 2 + 2^ cj f s + u; 2 

This transfer function is called the Disturbance Transfer. Figure 4-32 shows th' 
Bode plots of the disturbance transfer for two modal frequencies, one larger than 
the bandwidth of the filter and the other smaller. Note that for there 
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is actually an amplification of the unmodeled mode noise within the bandwidth of 
the filter, while for u/ f - < Wp there is an attenuation. 



Figure 4-32: Frequency Response of the Disturbance Transfer. 

There are four structural modes modeled into Detection Filter #3, two rigid 
modes and two bending modes. The rigid modes, with modal frequencies of 0 
rad/sec., totally reject the unmodeled mode disturbance at low frequencies, because 
the low frequency gain of the disturbance transfer is zero. On the other hand, the 
bending modes with modal frequencies of 11.4 and 31.4 rad/sec. have disturbance 
transfers with low frequency gains of 1.3 and 9.8 respectively. Therefore, the large 
error residuals in Fig. 4-30 are due to the bending modes. Specifically, the presence 
of the unmodeled modes generates neglectable low frequency errors in the residuals 
c* and corresponding to the rigid modes, significant low frequency errors in the 
residual £3 corresponding to the first bending mode and very large low frequency 
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errors in the residual cj” corresponding to the second bending mode. 

The output error vector € is calculated from e + as e — M e + . Therefore, the 
large residuals « 3 and produced by the unmodeled modes generate large 
unidirectional errors in the directions of c 3 and c 4 , respectively, in the output space 
of e. Hence, by choosing a different output transformation R, the errors in the 
directions c 3 and c 4 can be separated in such a way that they no longer affect the 
failure signature. The improved output transformation that achieves this 
separation is: 

R=(CAb 3 CAb 4 c 3 cj- 1 (4.28) 

and the new transformed output error vector E is: 


E = R € = 



where E 3 and contain the failure signatures of actuators 3 and 4 respectively, 
and E bl and E b2 only contain the error residuals due to the unmodeled modes, and 
correspond to the first and second bending mode respectively (actually E bl = t 
and E b2 = €+). 

Figure 4-33 shows the frequency response functions for the transfer from the 
unmodeled eighth mode to the residuals E 3 and E 4 obtained using the improved 
transformation R in (4.28), for Detection Filter #3 with no output filtering. These 


+ 00 
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response functions show a remarkable improvement when compared with the 
response functions in Fig. 4-28 obtained using the pseudo inverse as the output 
transformation. Note that the low frequency gains are actually zero. It is 
important to stress that these gains are zero because the low frequency gains of the 
disturbance transfers of the rigid modes are zero. That is, the zero low frequency 
gains in Fig. 4-33 are not particular to the uumodeled eighth mode, but they will 
occur in the transfers of all the unmodeled modes. 

The use of the improved output transformation solves the problem of the 
unmodeled modes only at low frequencies within the bandwidth of the filter. 
Output filtering is still needed to eliminate the high frequency components of the 
residuals that lie outside the bandwidth of the filter. Detection Filter#3 with 
output filtering and with the improved transformation was then tested in the 
presence of four unmodeled modes and failure of actuator 4 at T = 1 sec.. Figure 
4-34 shows the result of this simulation which was done under exactly the same 
conditions as the simulation in Fig. 4-30, but with a different output 
transformation. The improvement is remarkable. The low frequency components of 
the residuals have been totally eliminated and the failure signature is now clearly 
visible. The small error residuals consist mainly of a high frequency component at 
the frequency of the first unmodeled mode (i.e., 61.25 rad/sec.), that was not 
filtered out completely by the output filter. This high frequency component can be 
reduced further by increasing the order of the output filter. 

As mentioned earlier in this chapter, the state and output error equations 
with parameter errors in matrix C have the same structure as the error equations 
with unmodeled modes (see Fig. 4-11). Then, it is easy to see that the use of the 
improved output transformation produces transfers from to the residuals E 3 and 
E 4 with zero low frequency gains. Therefore, the use of the improved output 
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Figure 4-33: Detection Filter #3. Unmodeled Eighth Mode Transfer. 

R = Improved Transformation. 
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Figure 4-34: Detection Filter #3. Four Unmodeled Modes. Actuator 4 Failure 
at T = 1 sec.. R = Improved Transformation. Output Filtering. 

transformation also produces remarkable improvements in the performance of the 
filter in the presence of parameter errors in C. 

The parameter errors in the system matrix A consist only of the errors in the 
frequencies of the two bending modes. From Fig. 4-9 and Fig. 4-21, it can be seen 
that these errors affect only the loops corresponding to the bending modes, and 
therefore, only generate errors in the residuals £3 and Then, for the same 
reasons explained before, the improved transformation totally separates these errors 
from the failure signature. That is, the transfers from V, to E 3 and E 4 are zero at 
all frequencies, and therefore, the errors in the frequency of the two bending modes 
generate absolutely no errors in the residuals E 3 and E^. 

In the previous section, it was shown that Detection Filter #3 represented a 
large improvement in the performance in the presence of parameter errors in 
B. The new output transformation, however, only introduces a minor improvement. 
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With this transformation, only the errors in B corresponding to the modal shapes of 
the rigid modes generate errors in the residuals E 3 and E 4 . It is important to note 
that the shapes of the rigid modes are usually very well known. Therefore, the 
results shown in Figures 4-35 and 4-36 with 5% and 10% parameter errors in B are 
too pessimistic. 

Detection Filter #3, with output filtering and with the improved output 
transformation, was then tested in the presence of all types of model mismatch and 
with actuator 4 failure time at T = 1 sec.. Figure 4-35 shows the results of the 
simulation with four unmodeled modes and with 5% parameter uncertainties in all 
the matrices A, B and C. Figure 4-36 shows the results of the same simulation but 
with 10% parameter uncertainties in A, B and C. In both cases the failure 
signature is clearly visible, and therefoie, for the first time in this thesis it is 
possible to detect and identify the failure of actuator 3 and 4 under realistic 
conditions. 

The use of the new output transformation results in very large improvements 
in the performance of the filter in the presence of model mismatch. For this 
reason, in the rest of this section, the new output transformation will be analyzed a 
little further, and a new way of reading the detection information from the output 
error will be introduced. 

During nominal operation the error residuals due to model mismatch lie 
mainly in a plane spanned by the vectors c 3 and c 4 , as was shown earlier in this 
section. Then, when failure of actuator 3 or 4 occurs, the error residuals lie in the 
space spanned by the vectors CAb 3 , c 3 and c 4 or CAb 4 , c 3 and c 4 respectively. So, 
there are actually three output error spaces associated with this detection problem, 
one space corresponding to the nominal operation and two spaces corresponding to 
the failure of actuators 3 and 4. The improved output transformation produces 
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Figure 4-35: Detection Filter #3. Four Unmodeled Modes. Parameter Errors 
in A, B aDd C = b%. Actuator 4 Failure at T = 1 sec.. 

R = Improved Transformation. Output Filtering. 


excellent results because it is able to distinguish between these three spaces. 
During nominat operation, it generates large error residuals in E b ^ and E b . and 
when actuator 3 or 4 fails, it generates large error residuals in E 3 , E b ^ and E b ^ or 
E 4 , E b and E b respectively. Therefore, the residuals E 3 and E 4 can detect and 
identify the failure of actuator 3 and 4 in the presence of model mismatch. 


Detection Filter #3, as mentioned earlier in this section, produces 
unidirectional output with any actuator failure event. However, the improved 
transformation in (4.28) cannot be configured to detect more than two actuator 
failure events since two directions of the transformation have to be dedicated to 
suppress the effects of model mismatch. Therefore, Detection Filter #3 requires 
two improved output transformations in order to detect and identify the four 
actuator failure events. For example, one transformation can be configured to 
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Figure 4-30: Detection Filter #3. Four Unmodeled Modes, Parameter Errors 
in A, B and C = 10*%. Actuator 4 Failure at T = 1 sec.. 

R = Improved Transformation. Output Filtering. 

detect failure of actuators 1 and 2, and the other transformation can be configured 
to detect failure of actuators 3 and 4. Then, when actuator 1 fails, the residual Ej 
will grow while E,, will remain zero in the first transformation, and the residuals E 3 
and E 4 in the second transformation will both grow. Thus, the growth of three 
residuals are needed to identify the failed actuator. This is a very undesirable 
result because it complicates the decision algorithms and also increases the 
possibilities of incorrect identifications. 

A new way of reading the output error information that can identify the 
failure of any number of actuators, will now be explained. Even though this new 
way has some interesting advantages when compared with the output 
transformation, the main reason it is introduced here is to help gain more insight 
into the detection filter problem. 
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The new approach, used in conjunction with output filtering, is illustrated in 
Fig. 4-37. Note that the output error £ + is used instead of £, the reasons for this 
becoming clear later. Each component of the vector £ + is filtered with a bandwidth 
of 10 rad/sec. to eliminate the high frequency components produced by model 
mismatch (the output filters have to be identical to preserve the unidirectionality 
of £ + when a failure occurs). The filtered output error will be called v. When the 
jth actuator fails, the output error v remains in the single direction given by the 
vector btin (4.22). 



Figure 4-37: New Output Error Processing 

The new approach, unlike the transformations used so far, treats the 
detection and identification problems separately. Detection is based on the norm of 
the output error vector v. During nominal operation, the detection residual ||v|| is 
zero, if there is no model mismatch, and when a failure occurs, it becomes 
significantly nonzero. Therefore, by looking at the detection residual ||v||, the 
failure of any actuator (actually any component) can be detected. To identify 
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which component has failed, p directionality residuals are generated that represent 
the angles from 0 to QQ degress between the output error vector v and the output 
directions b+for j = 1, 2 . . p, corresponding to each actuator failure. During 
nominal operation, these residuals have no significance, and therefore, are not 
computed. When the jth actuator fails, the output error v remains in the direction 
of vector bj, and therefore the directionality residual is zero while the other 
directionality residuals have a nonzero value between 0 and 90 degrees. Hence, the 
actuator that failed is the one corresponding to the directionality residual that was 
zero (or closest to zero in the presence of model mismatch). 

In an operational system, the directionality residuals should be computed only 
when the detection residual becomes larger than a certain tolerance indicating that 
a component has failed. In this thesis, however, since the interest is to analyze the 
character of the residuals without actually making any decisions, the directionality 
residuals will be computed starting at T = 1 sec., which is the time at which the 
failure occurs. 

As mentioned earlier in this section, the presence of unmodeled modes 
generates large low frequency errors in the residuals and corresponding to the 
bending modes while generating negligible errors in the residuals t* and and 

£* can be eliminated by only using and an< ^ the components of b+ 
corresponding to the rigid modes, bj^j and bj^, in the computation of the detection 
and directionality residuals. But, since and c* are no longer needed in the 
residual generation, there is no need for the detection filter to estimate the 
amplitudes and amplitude rates of the bending modes and therefore, the states of 
the detection filter corresponding to the bending modes can be eliminated. Hence, 
the detection filter in the configuration shown in Fig. 4-26 will have only the four 
states corresponding to the rigid modes. However, the measurement transformation 
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matrix M must remain the same (4X4 nonsingular matrix) and therefore, the 
mode shapes of the first two bending modes are still needed in the design of this 
filter. 

A Detection Filter designed in this way is characterized as having 4 modeled 
modes and 2 detection and identification modes. Such a filter, then, has 4 states: 
the number of detection and identification modes times two. A detection filter 
with these characteristics was designed and was named Detection Filter #4. 

Detection Filter #4: 

# of modeled modes = 4 

jj. of detection and identification modes = 2 (4 states) 

Poles: 4 detection space poles at -lOrad/sec. 

Detection Filter #4 uses the errors in the estimates of the rigid modes 
amplitudes produced by the failure of an actuator to detect and identify the failed 
actuator. 

All the detection filters will be tested with output filtering, as shown in Fig. 
4-37, even in the case of no model mismatch, for the rest of this section. 

Figures 4-38 and 4-39 show the results of the simulations of Detection Filter 
#4 with no model mismatch and with failure of actuators 1 and 4, respectively, at 
T = 1 sec.. 

In both cases, the actuator failures create clearly evident signatures in the 
detection residuals. Also, in both cases, the directionality residuals corresponding 
to the failed actuators remain at zero while the other directionality residuals 
become clearly greater than zero. 

However, these simulations show a big difference in the ” ability” of Detection 
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Figure 4-39: Detection Filter #4. No Model Mismatch. 
Actuator 4 Failure at T = 1 sec.. Output Filtering. 


-114- 


Filter #4 to ideutify correctly the actuator that failed in the presence of model 
mismatch. It is intuitively clear that the ability of the detection filter to identify 
correctly the failure of the jth actuator in the presence of model mismatch is 
related to the smallest of the angular distances between the output direction 
corresponding to actuator j failure and the output directions corresponding to the 
other actuator failures. If this angle is large, it will take a large error residual due 
to model mismatch to "rotate" the output error in order to result in an incorrect or 
undetermined identification. If this angle is small, the opposite happens. From now 
on, this angle will be referred to as the "phase margin for actuator j* The phase 
margin for actuator 1 is 12 degrees while for actuator 4 is 44 degrees. Therefore, 
the ability of Detection Filter #4 to identify correctly the failure of actuator 1 in 
the presence of model mismatch, is much smaller than its ability to identify 
correctly the failure of actuator 4. 

The physical interpretation of the phase margin is very clear. The phase 
margin for actuator 1 corresponds to the angular distance between the failure 
direction of actuator 1 and the failure direction of actuator 2, as can be seen in Fig. 
4-38. These two actuators are placed next to each other, with a separation of 2 
feet, at the end of the beam (Fig. 3-1). The control effects of each of the actuators 
upon the rigid mode corresponding to the displacement of the beam, are identical. 
Therefore, what differentiates the failure of actuator 1 from the failure of actuator 
2, is the difference in the control of these actuators upon the rigid mode 
corresponding to rotation. Since they are placed close together at the end of the 
beam, they both have very similar influence on the rotation of the beam. 
Therefore, since the effect upon the rotation and displacement of the beam due to a 
control being applied to actuator 1 is very similar to the effect due to the same 
control being applied to actuator 2, it is very difficult to distinguish between the 
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failures of those two actuators from the amplitudes of the rigid modes. 

On the other hand, since actuator 4 is 3.5 feet away from the closest 
actuator, actuator 3, and since this actuator is at the center of the beam, the 
difference in the control of these actuators upon the rotation and translation of the 
beam, is much larger. This accounts for the large phase margin for actuator 4. 

Detection Filter #4 was then tested in the presence of four unmodeled modes. 
Figure 4-40 shows the results corresponding to the failure of actuator 1 and Fig 
4-41 shows the results corresponding to the failure of actuator 4. In both cases, the 
failure signatures in the detection residuals were much larger than the noise due to 
model mismatch. The dominant component in the noise has the frequency of the 
first unmodeled mode that was not completely eliminated by the output filter, as 
was the case for Detection Filter #3 in Fig. 4-34. The analysis of the directionality 
residuals is a little more involved. Whenever the failure signature is smaller than 
the noise due to model mismatch, the output error is dominated by the noise, and 
therefore, it is no longer unidirectional thus producing the large rapid vertical 
oscillations. Hence, the directionality residuals are significant only when the failure 
signature, in the detection residual, is larger than the background noise, and even 
in this case it is not guaranteed that the directionality residuals will identify 
correctly the failed actuator. For example, in the case of failure of actuator 1, the 
failure can be correctly, but marginally, identified for a very short period of time, 
even though the failure signature is much larger than the noise during most of the 
failure interval. On the other hand, the failure of actuator 4 can be clearly 
identified whenever the failure signature becomes a little larger than the noise. 

The poor performance of Detection Filter #4 when trying to identify the 
failure of actuator 1, was not due to large model mismatch poise, but rather due to 
the small phase margin for actuator I. In fact, the noise was much smaller than 
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Figure 4-40: Detection Filter #4. Four Unmodeled Modes 
Actuator 1 Failure at T = 1 sec.. Output Filtering. 
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the failure signature, allowing a very clear detection of the failure. One way of 
solving this problem is to increase the orders of the output filters in order to reduce 
the noise even further. Another solution, is to increase the phase margin for 
actuator 1 by augmenting the number of detection and identification modes. 
However, keep in mind that the estimation of the higher frequency modes are more 
affected by noise, and that was the reason for eliminating them in the design of 
Detection Filter #4. Therefore, this approach will be successful only if the 
improvement resulting from the larger phase margin is greater than the increase in 
noise due to model mismatch. To test this solution, Detection Filter #5 was 
designed with 4 modeled modes and 3 detection and identification modes. 

Detection Filter #5: 

# of modeled modes = 4 

# of detection and identification modes = 3 (6 states) 

#Poles: 6 detection space poles at -10rad/sec.. 

Detection Filter #5 uses the two rigid modes and the first bending mode in 
the process of detecting and identifying an actuator failure. Figure 4-42 shows the 
results of testing Detection Filter #5 with no model mismatch and with failure of 
actuator 1 at T = 1 sec.. Note that by incorporating the bending mode, the phase 
margin for actuator 1 increased from 12 to 37 degrees. 

Detection Filter #5 was then tested with four unmodeled modes and failure 
of actuator 1 at T = 1 sec.. The results of this simulation are shown in Fig. 4-43. 
The failure signature is still clearly evident in the detection residual and the result 
of adding a bending mode to the filter only produced a very small increase in the 
low frequency components of the noise. The larger phase margin, combined with a 
very low increase in noise, resulted in directionality residuals that clearly identify 
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Figure 4-42: Detection Filter #5. No Unmodeled Modes. 
Actuator 1 Failure at T = 1 sec.. Output Filtering 
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the failure of actuator 1. Therefore, it can be concluded that the addition of the 
bending mode to the filter was sucessful in improving the identification capabilities 
of the filter. 

To clarify what has been explained so far, the second bending mode will be 
added to the filter. Then, the resulting filter is Detection Filter #3 and the only 
difference is in the processing of the output error. Figure 4-44 shows the results of 
the simulation of Detection Filter #3 with no model mismatch and actuator 1 
failure at T = 1 sec.. The phase margin increased to 70 degrees, twice the phase 
margin of Detection Filter #5. Figure 4-45 shows the same simulation but with 
four unmodeled modes. The detection capabilities of the filter have been greatly 
impaired by the presence of large low frequency error residuals due to model 
mismatch. As explained earlier in this section, the large low frequency residuals are 
due to the large low frequency gain of the disturbance transfer corresponding to the 
second bending mode. However, note that the detection residual can still be 
generated from the rigid modes only, as was done in Detection Filter #4, while 
using the four modes in the calculation of the directionality residuals. In such a 
case, the advantage of the very low noise of the rigid modes in the detection of the 
failures, is combined with the advantage of the larger phase margin of the higher 
order filter in the identification of the failure. The results of the directionality 
residuals may be a little misleading. These residuals correctly indicate the failure 
of actuator 1. However, since the noise due to model mismatch is very slowly 
varying, it is not clear whether the residuals would continue to indicate correctly 
the failure of actuator 1 if the simulation was continued beyond 2 seconds. 
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Figure 4-43: Detection Filter #5. Four Uumodeled Modes. 
Actuator I Failure at T — 1 rec.. Output Filtering. 
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Figure 4-44: Detection Filter #3. No Model Mismatch. 
Actuator ! Failure at T = 1 sec.. Output Filtering. 
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Figure 4-45: Detection Filter #3. Four Unmodeled Moq q s. 
Actuator 1 Failure at T = 1 sec.. Output Filtering. 
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4.3.3 Detection Filter With Observer 

In the previous subsection, the problem of the large numbers of unmodeled 
nodes in the detection filter with N — p was solved by introducing a new output 
transformation. The large improvement that resulted from this was due to the 
rigid modes that rejected the low frequency noise, and therefore, provided clear 
detection and identifu ition information. In some cases, however, the model of the 
flexible structure might contain only one rigid mode or no rigid modes at all. In 
such a case, the noise rejection properties of the low frequency modes (see Eq. 
(4.27)) might nc, be good enough to reject the low frequency noise due to the 
unmodeled modes (the high frequency noise can always be eliminated with output 
filtering). The Detection Filter with Observer introduced in this subsection is 
a possible solution when this problem arises. However, it should be made clear that 
due to a lack of time the properties of such a configuration have not been analyzed 
in detail, and mainly the motivation for this approach is presented in this thesis. 

The problem of model mismatch in detection filters is that the requirement of 
having small error residuals seems to be incompatible with the requirement of 
having large failure signatures. The reduction in the error residuals due to model 
mismatch is achieved by improving the ability of the filter to track the output 
vector y, but in so doing, the ability of the filter to track the output vector y in 
the presence of an actuator failure will also improve, resuiting in a smaller failure 
signature. For example, the error residuals can be reduced by moving the detection 
space poles to the left in the s plane, but this also results in a smaller failure 
signature. The Detection Filter with Observer in Fig. 4-46 eliminates the error 
residuals due to those modes whose shapes and frequencies are known but are not 
modeled in the reduced order detection filter with N == p, without reducing the 
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Figure 4-40: The Detection Filter with Observer 
failure signature. 

The model of the observer contains all the structural modes of the beam that 
are known, while the model of the detection filter contains as many modes as the 
number of independent displacement measurements. The equations of the observer 
are 

i o = A o*o + B o\ + D o (y-y) 
y = C x„ 

and the equations of the detection filter are 


Xp — Ap Xp + Bp u c + Dp (yy - yp) 
yp = Cp Xp 
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where 


x 


o 



C 0 = [C F CJ 


The unmodeled modes of the filter - observer are defined as the modes that are not 
modeled either in the observer nor in the detection filter. The unmodeled modes of 
the detection filter are defined as the modes that are modeled in the observer but 
not modeled in the filter. The purpose of the observer is to obtain good estimates 
of the amplitudes of the unmodeled modes of the filter in the presence of model 
mismatch and more important, in the presence of actuator failures. These 
estimates are multiplied by their corresponding vectors c t from matrix C , and the 
results are subtracted from y to obtain yy. The detection filter, then, uses y F as 
the measurement vector. 

Therefore, when there is no model mismatch between the model of the beam 
and the model used in the observer and when there is no actuator failure, the 
estimates of the filter unmodeled modes are correct and the vector y F is the one 
resulting from those structural modes that are modeled in the filter. Hence, since 
the detection filter does not see any unmodeled modes, the residuals c are zero. 

When an actuator failure occurs, the result will depend on the ability of the 
observer to estimate the amplitudes of the unmodeled modes of the filter in the 
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presence of the actuator failure. If this ability is good, then the filter measurement 
vector y F will contain very little components resulting from the filter unmodeled 
modes, and therefore, the error residuals due to the filter unmodeled modes will be 
very small and identification of the failed component will be possible. On the other 
hand, if the observer is not good at estimating the amplitudes of the filter 
unmodeled modes in the presence of actuator failures, when an actuator failure 
occurs, y F will contain large components due to the errors in the estimates of the 
filter unmodeled mode amplitudes resulting in large error residuals in e that will 
make the identification impossible. 

Therefore, while the detection capabilities of the reduced order filter will be 
greatly improved by the use of any good observer, the identification capabilities 
will be improved only if this observer is good at estimating the amplitudes of the 
filter unmodeled modes in the presence of actuator failures. Note that this 
requirement is exactly the opposite of the requirement of the detection filter in 
which we want an actuator failure to produce large errors in the estimates. 
Furthermore, it seems that in this c;ise, there is no contradiction between requiring 
the observer to be insensitive to actvator failures and requiring the observer to be 
insensitive to model mismatch. 

A thorough analysis on whether it is possible to design obervers with such 
properties and how to design them has not been done in this thesis. It is 
reasonable to expect that observers that rely heavily on the measurements to 
produce the state estimates will have little sensitivity to actuator failures. An 
observer with this property has large bandwidth and high frequency poles. This is 
the case if a detection filter is used as the observer, since the failure signature 
caused by an actuator failure can be arbitrarily reduced by moving the detection 
space poles to the left in the s - plane. 
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The observer used in this thesis was designed using the continuous Kalman 
filter approach and a seven mode model of the beam. The resulting filter had some 
high frequency poles very close to the high frequency undamped structural modes, 
and therefore, it was necessary to use output filtering in the Detection Filter. 
However, it is believed that the output filtering can be avoided by a more careful 
design of the observer. 

Figure 4-47 shows the result of simulating the filter - observer combination 
with no model mismatch. That is, the evaluation model of the beam has seven 
modes, the observer has seven modes anl the reduced order detection filter has 
four modes. The reduced order detection filter used is Detection Filter #3 with the 
pseudo inverse as the output transformation. The initial conditions of the beam, 
observer and detection fitler are all zero and actuator 4 fails at T = 1 sec.. Figure 
4-47 should be compared with Fig. 4-30 (actually in Fig. 4-30 there is one more 
unmodeled mode but its contribution to the error residuals is negligible). The 
failure signature is now clearly visible. The use of the observer produced a 
remarkable improvement in the detection of the failure, as was expected, and also 
produced a very large improvement in the identification of the failure. 

The eighth structural mode was added to the evaluation model of the beam 
and the simulation repeated. The effect of the filter - observer unmodeled mode 
was so small that the output residuals were almost identical to the ones in Fig. 
4-47. 

Therefore, it can be concluded that the use of the observer was successful in 
improving the performance of Detection Filter #3 in the presence of unmodeled 
modes. The performance of the filter - observer configuration was not tested in the 
presence of parameter errors. 
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Figure 4-47: Detection Filter #3 with Observer. 
Failure of Actuator 4 at T = 1 sec.. Output Filtering. 
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Chapter 5 

Conclusions and Recommendations 

When the number of modes modeled into the detection filter was larger than 
the number of sensors, the detection filter could not be configured to detect more 
than three actuator failure events, and the corresponding state space was divided 
into detection and completion spaces. Moreover, the specification of the location of 
the filter poles did not result in a unique detector gain D. However, the algorithms 
used in this thesis, suggested by Beard, produced a unique gain matrix D by 
making the completion spaces invariant. As a result of all this, the filter design 
was complicated and relied heavily on software. 

The dynamics of these filters were very difficult to analyze and produced very 
unexpected results. The transfer properties of the detection filter that describe its 
performance in the presence of model mismatch, showed a surprising sensitivity to 
the location of the completion space poles explaining the results obtained in the 
simulations. When all of the filter poles were at -10rad/sec., the presence of an 
unmodeled mode created large low frequency error residuals that completely 
obscured the failure signature, and the presence of unrealistically low levels of 
parameter errors also resulted in unsuccessful detections and identification of the 
failure. By moving the completion space poles to -100rad/sec., the performance of 
the detection filter in the presence of model mismatch was remarkably improved. 
However, the detection filter still could not tolerate realistic levels of parameter 
errors. 

The detection filters with N = p (number of modeled structural modes equal 
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to number of sensors) have several advantages when compared to detection filters 
with N > p. Perhaps the most important of them is that detection filters with N 
= p can be configured to detect any number of actuator failure events by placing 
the filter poles in an appropriate way (e.g., placing all the poles at the same 
location.) The design of these filters is trivial and does not need any specialized 
software. Also, their dynamics are very simple and easy to analyze, since they do 
not involve completion spaces. This in turn allows the designer to make a thorough 
frequency response analysis of the effects of the different types of model mismatch, 
and make decisions on how to use the degrees of freedom that he has at his disposal 
like the location of the detection filter poles, bandwidth of the output filter, choice 
of output transformation, etc. Finally, another important advantage of these filters 
is that they perform much better in the presence of parameter errors. 

The problem of the detection filters with N = p is their large number of 
unmodeled modes. The high frequency components of the error residuals due to 
unmodeled modes lying outside the bandwidth of the filter can always be filtered 
out without reducing the failure signature. Then, the problem is caused by the 
components of the residuals lying within the bandwidth of the filter. This problem 
was dealt with in two different ways. First, the performance of the detection filter 
in the presence of model mismatch was improved by choosing a different output 
transformation. This resulted in a detection filter that was able to perform 
satisfactorily in the presence of the unmodeled modes and in the presence of 
realistic levels of parameter errors. The second approach was to use an observer to 
subtract the components of the unmodeled modes from the measurements y. 

From the analysis and simulations done in this thesis, it appears that a 
combination of all the techniques explained in subsection 4.3.1 has the best 
potential to solve the most severe problems of model mismatch in the detection of 
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actuator failures in flexible structures with displacement sensors. Such a 
combination refers to a detection filter with N = p, improved output 
transformation, output filtering and observer. 

The recommended areas of research related to the problem analyzed in this 
thesis are: 

1) to investigate whether different methods of assigning the n - completion 
poles, that do not require the invariance of the completion spaces, could not result 
in more robust detection filters, 

2) to investigate and develop different ways of reading the detection and 
identification information from the output error in the presence of model mismatch, 

3) to study all the aspects related to the detection filter with observer. For 
example, to study the design and properties of the observer that best fit this 
application, and to study the robustness properties of the filter - observer 
combination, 

4) to study the robustness properties of detection filters designed to detect 
actuator failures in flexible structures with both displacement and rate 
measurements, 

5) to study the performance of the detection filter in a more realistic 
environment that includes a control system and disturbances. 
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